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ABSTRACT 

We have performed a comprehensive chemical abundance analysis of the extremely metal -poor ([Ar/H]<- 
2) halo planetary nebula (PN) BoBn 1 based on IUE archive data, Subaru/HDS spectra, VLT/UVES archive 
data, and Spitzer/IRS spectra. We have detected over 600 lines in total and calculated ionic and elemental 
abundances of 13 elements using detected optical recombination lines (ORLs) and collisionally excited lines 
(CELs). The estimations of C, N, O, and Ne abundances from the ORLs and Rr, Xe, and Ba from the CELs 
are done the first for this nebula, empirically and theoretically. The C, N, O, and Ne abundances from ORLs 
are systematically larger than those from CELs. The abundance discrepancies apart from O could be explained 
by a temperature fluctuation model, and that of O might be by a hydrogen deficient cold component model. 
We have detected 5 fluorine and several slow neutron capture elements (the s-process). The amounts of [F/H], 
[Kr/H], and [Xe/H] suggest that BoBn 1 is the most F-rich among F detected PNe and is a heavy s-process 
element rich PN. We have confirmed dust in the nebula that is composed of amorphous carbon and PAHs with 
a total mass of 5.8 x 10~ 6 Mq. The photo-ionization models built with non-LTE theoretical stellar atmospheres 
indicate that the progenitor was a 1-1.5 Mq star that would evolve into a white dwarf with an ^0.62 Mq core 
mass and ^0.09 Mq ionized nebula. We have measured a heliocentric radial velocity of +191.6 ±1.3 km s -1 
and expansion velocity 2V exp of 40.5 ± 3.3 km s" 1 from an average over 300 lines. The derived elemental 
abundances have been reviewed from the standpoint of theoretical nucleosynthesis models. It is likely that the 
elemental abundances except N could be explained either by a 1.5 Mq single star model or by a binary model 
composed of 0.75 Mq + 1.5 Mq stars. Careful examination implies that BoBn 1 has evolved from a 0.75 Mq 
+ 1.5 Mq binary and experienced coalescence during the evolution to become a visible PN, similar to the other 
extremely metal -poor halo PN, K 648 in M 15. 

Subject headings: ISM: planetary nebulae: individual (BoBn 1, K 648), ISM: abundances, ISM: dust, Stars: 
Population II 



1. INTRODUCTION 

Planetary Nebulae (PNe) represent a stage in the evolution 
of low- to intermediate-mass stars with initial masses of 1-8 
Mq. At the end of their life, a star of such mass evolves first 
into a red giant branch (RGB) star, then an asymptotic giant 
branch (AGB) star, next a PN, and finally a white dwarf. Dur- 
ing their evolution, such stars eject a large amount of their 
mass. The investigation of chemical abundances in PNe en- 
ables the determination of how much of a progenitor's mass 
becomes a PN, when and how elements synthesized in the 
progenitor were brought to the surface, and how chemically 
rich the Galaxy was when the progenitors were born. 

Currently, over 1,000 objects are regarded as PNe in the 
Galaxy (Acker et al. 1991). Of these, about 14 objects have 
been identified as halo members from their location and kine- 
matics since the PN K 648 was discovered in M 15 (Pease 
1928). Halo PNe are interesting objects as they provide direct 
insight into the final evolution of old, low-mass halo stars, 
and they are able to convey important information for the 
study of low-mass star evolution and the early chemical con- 
ditions of the Galaxy. However, in extremely metal-poor and 

1 Space Telescope Science Institute, 3700 San Martin Dr., Baltimore, MD 
21218, USA; otsuka@stsci.edu 

2 Subaru Telescope, NAOJ, 650 North A'ohoku Place, Hilo, HI 96720, 
USA 

3 School of Science Education (Astronomy), Chungbuk National Univer- 
sity, 12 Gaeshin-dong Heungduk-gu, CheongJu, 

Chungbuk 361-763, Korea 

4 Okayama Astrophysical Observatory, NAOJ, 3037-5, Honjo, Kamogata- 
cho, Asakuchi-shi, Okayama, 719-0232, Japan 



C- and N-rich ([C,N/O]>0, [Ar/H]<-2) halo PNe, there are 
unresolved issues on chemical abundances and evolution time 
scales. BoBn 1 (PN G108.4-76.1) is one of the C- and N-rich 
and extremely metal-poor halo PNe ([C, N/0]>1, [Ar/H]=- 
2.22±0.09, [Fe/H]=-2.39±0.14; this work), which composes 
a class of PN together with K 648 (Otsuka 2007, see Table 
20) and H 4-1 (Otsuka et al. 2003). 

The progenitors of halo PNe are generally thought to be 
~0.8M Q stars, which is the typical mass of a halo star. Above 
mentioned three metal -poor C- and N-rich halo PNe, however, 
show signatures that they have evolved from massive progen- 
itors. For example, they would become N-rich, but would not 
C-rich if they have evolved from ~0.8 Mq single stars with 
[Fe/H]~-2.3 (Z=10~ 4 ), according to the current stellar evolu- 
tion models (e.g., Fujimoto et al. 2000). To become C-rich 
PNe, the third dredge-up (TDU) must take place in the late 
AGB phase. The efficiency of the TDU depends on the initial 
mass and composition, with increasing efficiency in models 
of increasing mass, or decreasing metallicity. At halo metal- 
licities, it is predicted that the TDU is efficient in stars with 
initial masses greater than ~1 Mq (Karakas 2010; Stancliffe 
2010). Also, current stellar evolutionary models predict that 
the post- AGB evolution of a star with an initial mass ~0.8M Q 
proceeds too slowly for a visible PN to be formed. The origin 
and evolution of halo PNe are still one of the unresolved big 
problems in this research field. 

How did these progenitor stars become visible C- and N- 
rich halo PNe? To answer this key question would deepen 
understanding of low-mass star evolution, in particular, ex- 
tremely metal-poor C-rich stars found in the Galactic halo, 
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and Galactic chemical evolution at early phases. If we can ac- 
curately estimate elemental abundances and ejected masses, 
then we can directly estimate elemental yields synthesized 
by PNe progenitors which might provide a constraint to the 
growth-rate of core mass, the number of thermal pulses and 
dredge-up mass. Hence we can build realistic stellar evolution 
models and Galactic chemical evolution models. We have ob- 
served these extremely metal-poor C- and N-rich halo PNe 
using the Subaru/High-Dispersion Spectrograph (HDS) and 
we also utilized collecting archival data carefully in order to 
revise our picture of these objects. In this paper, we focus on 
BoBn 1. 

The known nebular and stellar parameters of BoBn 1 are 
listed in Table 1 . Zijlstra et al. (2006) have associated BoBn 1 
with the leading tail of the Sagittarius (Sgr) Dwarf Spheroidal 
Galaxy, which traces several halo globular clusters. The helio- 
centric distance to the Sgr dwarf Galaxy is ^24.8 kpc (Kun- 
der & Chaboyer 2009), while the distance to this object is 
between 16.5 (Henry et al. 2004) and 29 kpc (Kingsburgh & 
Barlow 1992). 

BoBn 1 is a unique PN in that it might possess informa- 
tion about the chemical building-up history of the Galactic 
halo. Otsuka et al. (2008a) found that the [C/Fe] and [N/Fe] 
abundances of BoBn 1 are compatible with those of carbon- 
enhanced metal-poor (CEMP) stars. The C and N overabun- 
dances of CEMP can be explained by theoretical binary inter- 
action models (e.g., Komiya et al. 2007; Lau et al. 2007). 
Otsuka et al. (2008a) detected two fluorine (F) lines and 
found that BoBn 1 is the most F-enhanced and metal-poor PN 
among F-detected PNe. They found that the C, N, and F over- 
abundances of BoBnl are comparable to those of the CEMP 
star HE1305+0132 (Schuler et al. 2007). Through a com- 
parison between the observed enhancements of C, N, and F 
with the theoretical binary nucleosynthesis model by Lugaro 
et al. (2008), they concluded that BoBn 1 might share its ori- 
gin and evolution with CEMP-s stars such as HE1305+0132, 
and if that is the case the slow neutron capture process (the 
i-process) should be considered. 

According to current evolutionary models of low- to 
intermediate-mass stars, the ^-process elements are synthe- 
sized by slowly capturing neutrons during the thermal pulse 
AGB phase. The s-process elements together with carbon are 
brought to the stellar surface by the TDU. If we could find sig- 
natures that BoBn 1 has experienced binary evolutions such as 
mass transfer from a massive companion and coalescence, the 
issues on chemical abundances and evolutionary time scale 
would be simultaneously resolved. It would also be of great 
significance to reveal the origin of these elements in the early 
Galaxy through the study of metal-poor objects such as BoBn 
1 . We will search s-process elements and investigate their en- 
hancement in BoBn 1 . 

In this paper, we present a chemical abundance analysis of 
BoBn 1 using the newly obtained Subaru/HDS spectra, ESO 
VLT/UVES, Spitzer/IRS and WE archive data. We detect sev- 
eral candidate collisional excited lines (CELs) of s-process el- 
ements and optical recombination lines (ORLs) of N, O, and 
Ne. We determine ionic and chemical abundances of 13 ele- 
ments using ORLs and CELs. We construct a detailed photo- 
ionization model to derive the properties of the central star, 
ionized nebula, and dust. We also check consistency between 
our abundance estimations and the model. Finally, we com- 
pare the empirically derived abundances with the theoretical 
nucleosynthesis model values and discuss evolutionary sce- 
narios for BoBn 1. 



2. DATA & REDUCTIONS 
2. 1 . Subaru/HDS observations 

The spectra of BoBn 1 were taken using the High- 
Dispersion Spectrograph (HDS; Noguchi et al. 2002) attached 
to one of the two Nasmyth foci of the 8.2-m Subaru telescope 
atop Mauna Kea, Hawaii on October 6th 2008 (program ID: 
S08B-110, PI: M.Otsuka). In Fig. 1, we present the opti- 
cal image of BoBn 1 taken by the HDS silt viewer camera 
(~ 0."12 pixel" 1 , no filters) during the HDS observation. The 
sky condition was clear and stable, and the seeing was be- 
tween 0."4 and 0."6. The FWHM of the image is -1". BoBn 
1 shows a small protrusion toward the southeast. 

Spectra were taken for two wavelength ranges, 3600-5400 
A (hereafter, the blue region spectra) and 4600-7500 A (the 
red region spectra). An atmospheric dispersion corrector 
(ADC) was used to minimize the differential atmospheric dis- 
persion through the broad wavelength region. In these spec- 
tra, there are many recombination lines of hydrogen, helium, 
& metals and collisionally excited lines (CELs). These nu- 
merous spectral lines allowed us to derive reliable chemical 
compositions. We used a slit width of l."2 (0.6 mm) and a 
2x2 on-chip binning, which enabled us to achieve a nomi- 
nal spectral resolving power of /?=30000 with a 4.3 binned 
pixel sampling. The slit length was set to avoid overlap of the 
echelle diffraction orders at the shortest wavelength portion 
of the observing wavelength range in each setup. This corre- 
sponds to 8" (4.0mm), in which the nebula fits well and can 
allow us to directly subtract sky background from the object 
frames. The CCD sampling pitch along the slit length pro- 
jected on the sky is ~ 0."276 per binned pixel. The achieved 
S/N is >40 at the nebular continuum level even in both ends of 
each Echelle order. The resulting resolving power is around 
R >33 000, which derived from the mean of the full width 
at half maximum (FWHM) of narrow Th-Ar and night sky 
lines. All the data were taken as a series of 1800 sec expo- 
sure for weak emission-lines and 300 sec exposures for strong 
emission-lines. The total exposure times were 16200 sec for 
red region spectra and 7200 sec for blue region spectra. Dur- 
ing the observation, we took several bias, instrumental flat 
lamp, and Th-Ar comparison lamp frames, which were nec- 
essary for data reduction. For the flux calibration, blaze func- 
tion correction, and airmass correction, we observed a stan- 
dard star HR9087 at three different airmass. 

2.2. VLT/UVES archive data 

We also used archival high-dispersion spectra of BoBn 1, 
which are available from the European Southern Observatory 
(ESO) archive. These spectra were observed on August 2002 
(program ID: 069.D-0413, PI: M.Perinotto) and June 2007 
(program ID: 079.D-0788, PI: A.Zijlstra), using the Ultravio- 
let Visual Echelle Spectrograph (UVES; Dekker et al. 2000) 
at the Nasmyth B focus of KUEYEN, the second of the four 
8.2-m telescopes of the ESO Very Large Telescope (VLT) at 
Paranal, Chile. We call the August 2002 data "UVES1" and 
the June 2007 data "UVES2" hereafter. We used these data to 
compensate for unobserved spectral regions and order gaps in 
the HDS spectra. We normalized these data to the HDS spec- 
tra using the intensities of detected lines in the overlapped 
regions between HDS and UVES 1 & 2. 

These archive spectra covered the wavelength range of 
3300-6600 A in UVES 1 and 3300-9500 A in UVES2. The en- 
trance slit size in both observations was 1 1" in length and 1 ."5 
in width, giving R >30 000 derived from Th-Ar and sky lines. 
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TABLE 1 

Nebular and Stellar Parameters of BoBn 1 



Quantity 



Value 



References 



Name 

Position (J2000.0) 
Distance (kpc) 



Size (arcsec) 
logF(H/3) (erg cm" 2 s" 1 

com 

Rad. Velocity (km s" 1 ) 
Exp. Velocity (2V e x P ) 
T E (K) 
n € (cm -3 ) 
Abundances 



BoBn 1 (PNG108.4-76.1) 

a=0():37: 16.03 <5=-13:42:58.48 

22.5;29 

18.2;16.5 

24.8 

~2 (diameter) 
-12.54;-12.43 

-12.38;-12.53(observed);-12.44(de-redden) 
0.18;0.0;0.09 
191.6 (heliocentric) 
See Table 6 
See Table 8 
See Table 8 
See Table 19 



discovered by Boeshaar & Bond (1977) 

Hawley & Miller (1978);Kingsburgh & Barlow (1992) 

Mal'kov (1997);Henry et al. (2004) 

Kunder & Chaboyer (2009) 

This work (see Fig. 1) 

Cuisinier et al. (1996);Wright et al. (2005) 

Kwitter et al. (2003);This work;This work 

Cahn et al. (1992);Kwitter et al. (2003);This work 

This work 

This work 

This work 

This work 

therein Table 19 



log(L*/L ) 
log g (cm s~ 2 ) 
(M ) 

r*(K) 

Magnitude 



TABLE 2 

Journal of the HDS and UVES observations. 



Instr. Obs.Date seeing Range binning Exp. 

n (A) (sec) 

HDS 2008/10/06 0.4-0.6 3650-5400 2x2 1800x4 
0.4-0.6 3650-5400 2x2 600x3 
0.4-0.6 4600-7500 2x2 1800x9 
0.4-0.6 4600-7500 2 x 2 600 x3 
UVES 2002/08/04 0.8-1.5 3300-6600 lxl 2700x4 
2007/06/30 0.5-0.7 3300-9500 2x2 1500x2 

The CCDs used in UVES have 15 fim pixel sizes. For UVES 1 
an 1 x 1 binning CCD pattern was chosen. For UVES2 a 
2x2 on-chip binning pattern was chosen. The sampling pitch 
along the wavelength dispersion was ^0.015-0.02 A pixel" 1 
for UVES1 and ~0.03-0.04 A pixel" 1 for UVES2. The expo- 
sure time for UVES1 was 2700 sec x 4 frames, 10800 sec 
in total. The exposure time for UVES2 was 1500 sec x 2 
frames, 3000 sec in total. The standard star Feige 110 was 
observed for flux calibration. 

In Fig. 2, we present the combined HDS and UVES spec- 
trum of BoBn 1 normalized to the H/3 flux. The spectrum for 
the wavelength region of 3650-7500 A is from the HDS data 
and that of 3450-3650 A & >7500 A is from the UVES data. 

The observation logs are summarized in Table 2. The de- 
tected lines in the Subaru/HDS and VLT/UVES spectra are 
listed in Appendix A. 

2.3. IUE archive data 

We complemented optical spectra with UV spectra obtained 
by the International Ultraviolet Explorer (IUE) to derive C + , 
C 2+ , N 2+ , and N 3+ abundances from semi-forbidden lines 
C II], C III], N III], and N IV], since these emission lines can- 
not be observed in the optical region. These IUE spectra 
were retrieved from the Multi-mission Archive at the STScI 
(MAST). We collected the high- and low-resolution IUE spec- 
tra taken by the Short Wavelength Prime (SWP) and Long 
Wavelength Prime/Long Wavelength Redundant (LWP/LWR) 
cameras. Our used data set is listed in Table 3, and the 



Mal'kov (1997);Zijlstra et al. (2()06);This work 

Mal'kov (1997);This work 

Mal'kov (1997);This work 

Howard et al. (1997);Mal'kov (1997) 

This work 

Simbad data base 




2 0-8 



Off-Set (a) (arcsec) 

FIG. 1. — The optical image of BoBn 1 taken by the Subaru/HDS silt 
viewer camera. The seeing was 0.4" when we took this image and the nebular 
FWHM was measured to ~1". 

wavelength dispersion mode is indicated in the column 3. 
All the IUE observations were made using the large aperture 
(10.3x23 arcsec 2 ). SWP and LWP/LWR spectra cover the 
wavelength range of 1 150-1980 A and 1850-3350 A, respec- 
tively. For each SWP and LWP/LWR spectra, we did median 
combine to improve the S/N. The combined short wavelength 
spectrum was used to measure fluxes of emission-lines in 
< 19 10 A because this allowed us to separate C III] A 1906/08 
and C iv A1548/51 lines. C m]A1906/08 are important as 
a density diagnostic. The combined long wavelength spec- 
trum was for measurements of emission-line fluxes in >2000 
A. The measured line fluxes were normalized to the H/3 
flux using theoretical ratios of He II 7(A1640)/(A4686) for the 
short wavelength spectrum and /(A25 12)/(A4686) for the long 
wavelength spectrum, respectively, adopting an electron tem- 



3.57;3.72;3.07 
5.52;6.5 
0.575;0.62 
125 000;96300 
125260 

16(B),14.6(R),16.13(J),15.62(H),15.18(K) 
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FIG. 2.— The spectrum of BoBn 1. The flux density is normalized so that H/3 flux F(Hfi) = 100. 
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TABLE 3 
Journal of IUE observations. 



Camera Data ID. disp. 



Range 
(A) 



Obs.Date 



Exp.time 

(sec) 



LWR 


16515 


low 1850-3350 1983-08-03 


6780 


LWP 


23692 


low 1850-3350 1992-08-13 


1500 


LWP 


23697 


low 1850-3350 1992-08-14 


7200 


LWP 


23699 


low 1850-3350 1992-08-15 


12000 


SWP 


45367 


high 1150-1980 1992-08-18 


7200 


LWP 


23713 


low 1850-3350 1992-08-18 


1800 


SWP 


45369 


high 1150-1980 1992-08-18 


10500 


SWP 


45371 


high 1150-1980 1992-08-19 


19 800 


SWP 


45386 


high 1150-1980 1992-08-21 


9000 



perature T e = 8840 K and density n e = 10 cm as given by 
Storey & Hummer (1995), then normalized to the H/3 flux. 
The interstellar extinction correction was made using equa- 
tion (1) (see section 3.1). The observed and normalized fluxes 
of detected lines are listed in the columns 4 and 5 of Table 4, 
respectively. 



2.4. Spitzer archive data 

We used two data sets (program IDs: P30333 PI: A.Zijlstra; 
P30652 PI: J.Bernard-Salas) taken by the Spitzer space tele- 
scope in December 2006. The data were taken by the Infrared 
Spectrograph (IRS, Houck et al. 2004) with the SH (9.5-19.5 
lira), LH (5.4-37/zm), SL (5.2-14.5 H and LL (14-38/an) 
modules. In Fig. 3 we present the Spitzer spectra of BoBn 
1 . We downloaded these data using Leopard provided by the 
Spitzer Science Center. The one-dimensional spectra were 
extracted using Spice version cl5.0A. We extracted a region 
within ±1" from the center of each spectral order summed up 
along the spatial direction. For SH and LH spectra, we sub- 
tracted sky background using off-set spectra. We normalized 
the SL and LL data to the SH and LH using the measured 
fluxes of [S iv] A10.5^m, H I A12.4 ^m, [Ne n] A12.8/xm, 
[Ne III] A15.6^m, and [Ne III] A36.0^m. Finally, the mea- 
sured line fluxes were normalized to the H/3 flux. The ob- 
served line ratio H 1 1(1 1.2^m)/7(A4861) (3.1 x 10" 3 ) is con- 
sistent with the theoretical value (3. 15 x 10" 3 ) for T e = 8840 K 
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and n e = 10 4 cm -3 as given by Storey & Hummer (1995). We 
did not therefore perform interstellar extinction correction. 

The observed and normalized fluxes of detected lines are 
listed in Table 5. In addition to the ionized gas emissions, 
the amorphous carbon dust continuum and the polycyclic aro- 
matic hydrocarbons (PAHs) feature around 6.2, 7.7, 8.7, and 
1 1 .2 fim are found for the first time. In Fig. 4, we present 
these PAH features. The 11.2 fim emission line is a complex 
of the narrow width H I 11.2 /itm and the broad PAH 11.2 
fim. The 6.2, 7.7, and 8.7 jim bands emit strongly in ionized 
PAHs, while the 11.2 fim does in neutral PAHs (Bernard-Salas 
et al. 2009). According to the PAH line profile classifications 
by Peeters et al. (2002) and van Diedenhoven et al. (2004), 
BoBn l's PAH line-profiles belong to class B. Bernard-Salas 
et al. (2009) classified 10 of 14 Magellanic Clouds (MCs) 
PNe into class B based on Spitzer spectra. In measuring PAH 
band fluxes, we used local continuum subtracted spectrum by 
a spline function fitting. We followed Bernard-Salas et al. 
(2009) and measured integrated fluxes between 6.1 and 6.6 
fim for the 6.2 fim PAH band, 7.2-8.3 /jm for the 7.7 ^m PAH, 
8.3-8.9 fim for the 8.6 fim PAH, and 1 1 . 1 - 1 1 .7 fim for the 1 1 .2 
fim. The observed PAH flux ratios 7(6.2/im)//(l 1.2/im) and 
I{1 .1 fim)II{\ 1.2/im) follow a correlation among MCs PNe, 
shown in Fig. 2 of Bernard-Salas et al. (2009). BoBn 1 has a 
hot central star (>10 5 K), so that ionized PAH might be dom- 
inant. However, these line ratios of BoBn 1 are somewhat 
lower than those of excited MCs PNe. One must take a look 
at a part of the neutral PAH emissions in a photodissociation 
region (PDR), too. 

We found a plateau between 10 and 14 fim, which are 
believed to be related to PAH clusters (Bernard-Salas et al. 
2009). Meanwhile, MgS feature around 30 /im sometimes 
observed in C-rich PNe, was unseen in BoBn 1 . 

2.5. Data reduction 

Data reduction and emission line analysis were performed 
mainly with a long-slit reduction package noao . t wodspec 
in I RAF 5 . Data reduction was performed in a standard man- 
ner. 

First, we made a zero-intensity level correction to all frames 
including flat lamp, object, and Th-Ar comparison frames 
using the overscan region of each frame and the mean bias 
frames. We also removed cosmic ray events and hot pixels 
from the object frames. Second, we trimmed the overscan 
region and removed scattered light from the flat lamp and ob- 
ject frames. Third, we made a CCD sensitivity correction to 
the object frames using the median flat frames. Fourth, we ex- 
tracted a two-dimensional spectrum from each echelle diffrac- 
tion order of each object frame and made a wavelength cali- 
bration using at least two Th-Ar frames taken before and after 
the object frame. We referred to the Subaru/HDS comparison 
atlas 6 and a Th-Ar atlas. For the wavelength calibration, we 
fitted the wavelength dispersion against the pixel number with 
a fourth- or fifth-order polynomial function. With this order, 
any systematic trend did not show up in the residuals and the 
fitting appears to be satisfactory. We also made a distortion 
correction along the slit length direction using the mean Th- 
Ar spectrum as a reference. We fitted the slit image in the Th- 
Ar spectrum with a two-dimensional function. For the HDS 

5 I RAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in Astron- 
omy (AURA), Inc., under a cooperative agreement with the National Science 
Foundation. 

6 http://www.naoj.org/Observing/Instruments/HDS/wavecal.html 



TABLE 4 

The detected lines in the IUE spectra. 



A kb Ion 


/(A) 


F(A) 


/(A) 


(A) 




(erg s" 1 cnT 2 ) 


[7(H/3)=100] 



1485 N IV] 1.306 1.41(-13)±5.47(-14) 45.81 ± 17.84 

1548 Civ 1.239 3.27(-12) ± 5.04(-14) 1052.6 ± 20.22 

1551 Civ 1.237 1.62(-12) ± 4.27(-14) 519.3 ± 14.96 

1640 Hell 1.177 5.13(-13) ± 7.40(-14) 162.97 ± 23.57 

1750 N III] 1.154 1.52(-13) ±4.16(-14) 48.06±13.16 

1906 Cm] 1.255 2.60(-12) ± 2.50(-14) 838.88 ± 12.65 

1908 Cm] 1.258 1.87(-12) ± 2.33(-14) 602.72 ± 10.28 

2324 C II] 1.388 2.87(-13) ± 3.60(-14) 36.64 ± 4.63 
+[Om] 

2424 [Ne IV] 1.134 1.41(-13) ± 1.46(-14) 17.12 ±1.78 

2512 Hell 0.969 2.66(-14) ± 1.29(-14) 3.13 ± 1.52 

Note. — X(-Y) stands for X x 10" r . /(H/3) = 3.63(-13) ± 6.47(-14) 
(see section 3.1). 



TABLE 5 

The detected lines in the Spitzer spectra. 



Aiab. Ion 


F(A) 


/(A) 


(p,m) 


(erg s _1 cnT 2 ) 


[/(H/3)=100] 



6.2 
7.7 
8.6 
10.5 
11.3 
11.3 
12.4 
12.5 
12.8 
14.6 
15.6 
16.4 
18.7 
19.1 
25.9 
36.0 



PAH 

PAH 

PAH 

[SIV] 

Hi 

PAH 

Hi 

He I 

[Ne II] 

He I 

[Ne III] 

He I 

[Sill] 

Hi 

[Oiv] 
[Ne III] 



62(-14) ± 1.52(-15) 
13(-14) ±2.31(-15) 
02(-14) ± 1.23(-15) 
36(-15) ± 1.93(-16) 
180-15) ± 1.810-16) 
97(-14) ± 9.650-15) 
89(-15) ± 3.0K-16) 
72(-15) ±4.17(-16) 
53(-15) ±3.13(-16) 
85(-16) ± 2.40(-16) 
17(-13) ±4.98(-15) 
190-15) ± 2.65(-16) 
650-15) ± 1.790-16) 
140-15) ± 2.64(-16) 
77(-14) ± 5.20(-16) 
100-14) ± 7.04(^15) 



7.22 ± 1.35 
22.39 ± 4.04 
5.57 ± 1.05 
1.92 ±0.05 
0.31 ±0.05 
13.68 ±3.61 
1.02 ±0.08 
0.45 ±0.11 
2.49 ± 0.08 
0.15 ±0.06 
161.13 ± 1.30 
0.57 ± 0.07 
0.69 ± 0.05 
0.30 ± 0.07 
12.46 ±0.14 
13.32 ± 1.84 



Note. — X(- 
(see section 3.1). 



Y) stands for X X 10- F . /(H/3) = 3.63(-13) ± 6.47(-14) 



spectra, we adopted fourth- and third-order polynomial func- 
tions for the wavelength and space directions, respectively. 
The fitting residual was of the order of 10" 4 A. For the UVES 
spectra, we adopted third- and second-order polynomial func- 
tions for the wavelength and space directions, respectively. 
The fitting residual was of the order of 10" 3 A. Fifth, we de- 
termined a sensitivity function using sky-subtracted standard 
star frames and obtained sky-subtracted and flux-calibrated 
two-dimensional PN spectra. The probable error in the flux 
calibration was estimated to be less than 5 %. Finally, we 
made a spatially integrated one-dimensional spectrum, and 
we combined all the observed echelle orders using I RAF task 
scombine. 

In measuring emission line fluxes, we assumed that the line 
profiles were all Gaussian and we applied multiple Gaussian 
fitting techniques. 

3. RESULTS 

3.1. Interstellar reddening correction 

We have detected over 600 emission lines in total. Before 
proceeding to the chemical abundance analysis, it is neces- 
sary to correct the spectra for the effects of absorption due to 
the Earth's atmosphere and interstellar reddening. The for- 
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FIG. 3. — The Spitzer spectra of BoBn 1 obtained by the SH and LH (upper panel) and the SL and LL modules (lower panel). 
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mer was performed using experimental functions measured 
at the Keck observatories and the ESO/VLT. The interstellar 
reddening correction was made by determining the redden- 
ing coefficient at H/3, c(Hf3). We fitted the observed inten- 
sity ratio of Ha to H/3 with the theoretical ratios computed 
by Storey & Hummer (1995). Two different situations are 
assumed for some lines: Case A assumes that the nebula is 
transparent to the lines of all series of hydrogen; Case B as- 
sumes the nebula is partially opaque to the lines of the Lyman 
series but is transparent for the Balmer series of hydrogen. 
Initially, we assumed that T e = 10 4 K and n e = 10 4 cm" 3 in 
Case B, and we estimated c(H/3) = 0.087 ± 0.004 from the 
HDS spectra. That is an intermediate value between Cahn 
et al. (1992) and Kwitter et al. (2003) (see Table 1). For 



the UVES spectra, we estimated c(H/3) = 0.066 by the same 
manner. From Seaton's (1979) relation c(H/3) = lATE(B-V) 
one obtains E(B-V) = 0.06 for the HDS spectra and 0.04 for 
UVES spectra, which are comparable to the Galactic value 
(0.02) to the direction to BoBn 1 measured by the Galactic 
extinction model of Schlegel et al. (1998). 

All of the line intensities were then de-reddened using the 
formula: 



log 



10 



/(A) 



/(H/3) 



= log 



10 



F(X) 



F(H/3) 



+ c(H/3)/(A), (1) 



where /(A) is the de-reddened line flux; F(X) is the observed 
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line flux; and /(A) is the interstellar extinction at A. We 
adopted the reddening law of Cardelli et al. (1989) with the 
standard value of R v = 3.1 for /(A). We observed F(Hf3) = 
2.56xl0" 13 ± 2.13xl(T 16 erg s" 1 cm" 2 (X(-Y) stands for X 
x 10~ F , hereafter) within the 1 ."2 slit in the HDS observation. 
We estimated captured light from BoBn 1 using the image 
presented in Fig. 1, to be about 86.8 % of the light from BoBn 
1 in the HDS observation. The intrinsic observed H/3 flux is 
2.95(-13) ± 3.45(-16) erg s" 1 cm" 2 and the de-reddened U/3 
flux is 3.63(— 13) ± 6.47(-14) erg s -1 cm" 2 including the error 
ofc(H^). 

3.2. Radial and expansion velocities 

We present the line-profiles of selected ions in Fig. 5. The 
observed wavelength at the time of observation was corrected 
to the averaged line-of-sight heliocentric radial velocity of 
+ 191.60 ± 1.25 km s" 1 among over 300 lines detected in the 
HDS spectra. The line-profiles can be represented by a single 
Gaussian for weak forbidden lines such as [Ne v] A3426 and 
metal recombination lines such as O II A4642. For the oth- 
ers, the profiles can be represented by the sum of two or three 
Gaussian components. 

Most of the detected lines are asymmetric profile, in partic- 
ular the profiles of low-ionization potential ions show strong 
asymmetry. The asymmetric line-profiles are sometimes ob- 
served in bipolar PNe having an equatorial disk structure. The 
similar line-profiles are also observed in the halo PN H4-1 
(Otsuka et al. 2006). H4-1 has an equatorial disk structure and 
multi-polar nebulae. The elongated nebular shape of BoBn 1 
(Fig. 1) might indicate the presence of such an equatorial 
disk. The receding ionized gas (especially, low-ionization po- 
tential ions) from the observers around the central star would 
be strongly weakened by the equatorial disk. In contrast, the 
relatively large extent bipolar flows perpendicular to the equa- 
torial disk might be unaffected by the disk. Due to such a 
geometry, we observe asymmetric line-profiles. 

In Table 6, we present twice the expansion velocity 2V exp 
measured from selected lines. When we fit line-profile with 
two or three Gaussian components, we define that 2Ve Xp corre- 
sponds to the difference between the positions of the red and 
blue shifted Gaussian peak components. We call expansion 
velocity measured by this method '2V exp (a)'. 

When we can fit line-profile with single Gaussian, we de- 
termine 2Ve Xp from following equation, 

2 Vexp = (^FWHM ~~ ^therm ~~ ^instr) ^ ) (2) 

where Vfwhm is the velocity FWHM of each Gaussian compo- 
nent. Vtherm and Vi nst r (~9 km s" 1 ) are the thermal broadening 
and the instrumental broadening, respectively (Robinson et al. 
1982). Vtherm is represented by 2\A{T/[/A) x / 2 , where T4 is the 
electron temperature (in units of 10 4 K) and A is the atomic 
weight of the target ion. For CELs, we adopted T4 listed in 
Table 9. For ORLs, we adopted T 4 = 0.88. We call twice the 
expansion velocity measured by Equation (2) '2Ve Xp (b)\ We 
converted 2V exp (a) into 2V eX p(b) using the relation 2V exp (b) = 
(1.74 ± 0.12) x 2V exp (a) for 2V exp (b) > 16 km s" 1 , which 
can be applied only to BoBn 1. We adopted 2V exp (b) as twice 
the expansion velocities for BoBn 1 . The averaged 2V exp (b) 
is 40.5 ± 3.28 km s" 1 among selected lines listed in Table 6 
and 33.04 ± 2.61 km s" 1 among over 300 detected lines in the 
HDS spectra. 

In Fig. 6 we present relation between 2V exp (b) and ioniza- 
tion potential (I. P.). When we assume that BoBn 1 has a stan- 



TABLE 6 

Twice the expansion velocities from selected lines . 



Ion A lab LP. 2V exp (a) 2V exp (b) 
(A) (eV) (kms' 1 ) (fans' 1 ) 



[Ne V] 


3425 


97.1 




10.0 


Ne II 


3694 


41.0 




34.4 


[Oil] 


3726 


13.6 


32.0 


55.6 


[Ne III] 


3868 


41.0 




41.5 


[FlV] 


4UoU 


62. / 




1 J.U 


C III 


4187 


47.9 




31.9 


C II 


4267 


24.4 




48.4 


N III 


4379 


47.5 




40.9 


Nil 


4442 


29.6 




41.8 


11 


4642 


35.5 




24.2 


He II 


4686 


54.4 




26.9 


[Ariv] 


4711 


40.7 




32.3 


[Ne IV] 


4724 


63.5 




8.9 


[Fll] 


4798 


17.4 




33.3 


H/3 


4861 


13.5 




45.0 


[Fe III] 


4881 


16.2 




32.6 


rn ttti 
[U IIIJ 


DUU / 


jj.j 




/to c 


[Nl] 


5198 





35.3 


61.3 


[CI III] 


5517 


23.8 




56.3 


Civ 


5812 


64.5 




14.3 


[Fill] 


5733 


35.0 




48.6 


He 1 


5876 


24.6 




32.7 


[Ol] 


6300 





48.0 


83.3 


[Nil] 


6548 


14.5 


37.7 


65.4 


[SII] 


6716 


10.4 


32.1 


55.7 


[Fe IV] 


6741 


30.7 




53.4 


[Arm] 


7135 


27.6 




48.0 


[S ml 


9069 


23.3 




49.7 



NOTE. — The probable error of 2V ex p(b) is within 5 km s We assume 
2Vex P (b) as twice the expansion velocity of BoBn 1 (see text). 

dard ionized structure (i.e., high LP lines are emitted from 
close regions to the central star and low LP lines are from 
far regions), expansion velocity of BoBn 1 seems to be pro- 
portional to the distance from the central star. BoBn 1 might 
have Hubble type flows. We found that 2V exp (b) values from 
ORLs are slightly smaller than CELs with the same LP, for 
example, Oil & [Om] (35.5 eV) and Nell & [Nelll] (41.0 
eV). O II and Ne II might be emitted from colder regions than 
[O III] and [Ne III] do. 

3.3. Plasma diagnostics 

3.3.1. CEL diagnostics 

We have detected a large number of collisionally excited 
lines (CEL), useful for estimations of the temperatures (T e ) 
and densities (n e ). The electron temperature and density diag- 
nostic lines analyzed here arise from various ions, which have 
a wide variety of ionization potentials ranging from ([N I] 
& [O I]) to 63.5 eV ([Ne IV]). We have examined the electron 
temperature and density structure within the nebula of BoBn 
1 using 17 diagnostic line ratios. The [Ol], [Nelll], [Neiv], 
[S II], and [S III] zone electron temperatures and [N I], Cm] 
and [Ne III] zone electron densities are estimated for the first 
time. Electron temperatures and densities were derived from 
each diagnostic ratio for each line by solving level populations 
for a multi-level (^5 for almost all the ions) atomic model 
using the collision strengths (J (j>i) and spontaneous tran- 
sition probabilities Ajj for each ion from the references given 
in Table 7. 

The derived electron temperatures and densities are listed 
in Table 8. Fig. 7 is the diagnostic diagram that plots the loci 
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FIG. 5. — The line-profiles of selected ions. Vertical and horizontal axes are scaled flux density and velocity with respect to the systemic radial velocity of 
+191.60 km s" 1 , respectively. All are from the HDS spectra except [NeV]A3426 and [S III] A9069 which are from UVES1 and UVES2, respectively. 



of the observed diagnostic line ratios on the n e -T e plane. This 
diagram shows that most CELs in BoBn 1 are emitted from 
7W2000-16000 K and log 10 « e ~3.5 cm" 3 ionized gas. 

First, we calculated electron densities assuming a constant 
electron temperature of 12 800 K. Estimated electron densities 
range from 1030 ([Nl]) to 5740 cm" 3 ([S II]). Although [S II] 
and [O n] have similar ionization potentials, a large discrep- 
ancy between their electron densities is found (see Table 8). 
Kniazev et al. (2008) and Kwitter et al. (2003) estimated [S II] 



electron densities as large as 9600 and 7100 cm , respec- 
tively. Stanghellini & Kaler (1989), Copetti & Writzl (2002), 
and Wang et al. (2004) found that the [S II] density is sys- 
tematically larger than the [O II] density in a large number of 
samples. The curve yielded by the [S ll]A6716/31 ratio in the 
n e vs T e plane indicates higher electron density than critical 
density of these lines, 1600 & 4100 cm -3 at T e = 12 800 K for 
A6716 and A6731, respectively (cf. 5740 cm" 3 in Fig. 7). This 
density discrepancy is not due to the errors in the [O II] atomic 
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TABLE 7 

Atomic data references for CELs. 



80 



60 



£ 40 - 



w 20 



Line transition probabilities A 



collisional strength f2;. 



[Cl] 

[Cii] 

Cm] 

Civ 

[Nl] 

[Nil] 

NlII] 

NIV] 

[Ol] 

[On] 

[Oin] 

[Oiv] 

[FII] 

[Fill] 

[FIV] 

[Ne II] 

[Ne III] 

[Ne IV] 

[Ne V] 

[S II] 

[S in] 

[SIV] 

[Clm] 
[Cliv] 

[Arm] 

[ArlV] 

[Fe III] 

[FelV] 

[KrlV] 

[Krv] 

[RbV] 

[Xe III] 

Ban 



Froese-Fischer & Saha (1985) 

Nussbaumer & Storey (1981); Froese-Fischer (1994) 

Wieseetal. (1996) 

Wieseetal. (1996) 



(1996) 
(1996) 

(1995); Froese-Fischer (1983) 

(1996) 

(1996) 



Wiese et al. 
Wiese et al. 
Brage et al 
Wiese et al. 
Wiese et al. 
Wieseetal. (1996) 
Wieseetal. (1996) 
Wiese etal. (1996) 

Storey & Zeippen (2000); Baluja & Zeippen (1988) 
Naqvi(1951) 

Garstang (1951); Storey & Zeippen (2000) 

Saraph&Tully (1994) 

Mendoza (1983); Kaufman & Sugar (1986) 

Becker et al. (1989); Bhatia & Kastner (1988) 

Kaufman & Sugar (1986); Bhatia & Doschek (1993) 

Verner et al. (1996); Keenan et al. (1993) 

Tayal & Gupta (1999) 

Johnson et al. (1986); Dufton et al. (1982); Verner et al. 
Mendoza & Zippen (1982a); Kaufman & Sugar (1986) 
Mendoza & Zippen (1982b); Ellis & Martinson (1984) 
Kaufman & Sugar (1986) 
Mendoza (1983); Kaufman & Sugar (1986) 
Mendoza & Zippen (1982a); Kaufman & Sugar (1986) 
Garstang (1957); Nahar & Pradhan (1996) 
Froese-Fischer & Rubin (1998); Garstang (1958) 
Biemont & Hansen (1986) 
Biemont & Hansen (1986) 
Persson et al. (1984) 
Biemont et al. (1995) 

Klose et al. (2002) 



Johnson et al. (1987); Pequignot & Aldrovandi (1976) 
Blum & Pradhan (1992) 
Berrington et al. (1985) 

Badnell & Pindzola (2000); Martin et al. (1993) 

Pequignot & Aldrovandi (1976); Dopita et al. (1976) 

Lennon & Burke (1994) 

Blum & Pradhan (1992) 

Ramsbottom et al. (1994) 

Bhatia & Kastner (1995) 

McLaughlin & Bell (1993); Pradhan (1976) 

Lennon & Burke (1994) 

Blum & Pradhan (1992) 

Butler & Zeippen (1994) 

See Text 

Lennon & Burke (1994) 
Saraph & Tully (1994) 
McLaughlin & Bell(2000) 
Ramsbottom et al. (1998) 
Lennon & Burke (1994) 
Ramsbottom et al. (1996) 
Froese Fischer et al. (2006) 
(1996) Dufton et al. (1982) 

Ramsbottom et al. (2001) 
Galavis et al. (1995) 

Galavis et al. (1995) 
Zeippen et al. (1987) 
Zhang (1996) 
Zhang & Pradhan (1997) 
Schoning (1997) 
Schdning (1997) 

Schoning & Butler (1998) 
Schoning & Butler (1998) 
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FIG. 6. — Relation between 2V exp (b) and ionization potential (LP.). The 
filled circles are the values from the CELs and the open circles are from the 
ORLs. 

data. Wang et al. (2004) also found the density discrepancy 
between [S II] and [O II] that might be likely caused by er- 
rors in the transition probabilities of [O II] given by Wiese 
et al. (1996). In the case of BoBn 1, this possibility can be 
ruled out because we obtained similar [O II] electron densities 
even when with the other transition probabilities. The high 
[S II] density might be due to high-density blobs in the outer 
nebula. This could have contributed to producing the small 
[S II]A6717/A6731 ratio and give rise to an apparently high 
density. Zhang et al. (2005a) pointed out the possibility that 
a dynamical plow by the ionization front effects yields large 
density of [S II] because the ionization potential of S + is close 



to the H + edge. Since the estimated upper limit to the [S II] 
density is close to the H + density derived from the Balmer 
decrement (see below), this explanation might be plausible. 
In BoBn 1, caution is necessary when using the [S II] electron 
density. 

Next, we calculated the electron temperature. An average 
electron density of 3370 cm" 3 , which excluded the « £ ([Nl]) 
and n £ ([S II]), was adopted when estimating electron temper- 
atures except the r e ([Ol]). [Nl] and [Ol] are representative 
of the very outer part of the nebula and probably do not co- 
exist with most of the other ions. r e ([Ol]) was, therefore, 
estimated using n e ([N I]). 

To obtain the [Nil], [Oil] and [Olll] temperatures it is 
necessary to subtract the recombination contamination to the 
[N II] A5755, [O II] AA7320/30, and [O III] A4363 lines, re- 
spectively. For [N II] A5755, Liu et al. (2000) estimated the 
contamination to [Nil] A5755, I R ([N ll]A5755) in the range 



5000 ^ r e ^ 20 000 K as 



/ R ([NII]A5755) 



/(H/3) 



: 3.19 



lov 



0.33 



2+ 



H+ 



(3) 



where N 2+ /H + is the doubly ionized nitrogen abun- 
dance. Adopting the value derived from the ORL analy- 
sis (see Section 3.6) and using Equation (3), we estimated 
/ s ([Nll]A5755) -0.1, which is approximately 7 % of the ob- 
served value. Given the corrected [N II] A5755 intensity, the 
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[N II] temperature is 12 000 K, which is 400 K lower than that 
obtained without taking into account the recombination effect. 

The same effect also exists for the [O II] AA7320/30 lines. 
We estimated the recombination contribution using the dou- 
bly ionized oxygen abundance derived from O II lines and the 
equation of Liu et al. (2000) for these lines in the range 5000 
^ T e ^ 10000 K, 



7 R ([OII]AA7320/30) 



1± 
io 4 



0.44 



o 



2+ 



9-36 ( — I x — (4) 



Using Equation (4), we estimated a contribution of ^7 % of 
the observed value and obtained 12 100 K, lower by 700 K 
than that without the recombination contribution. For [O III] 
A4363, we estimated the recombination contribution using 
the 3+ abundance derived from the fine-structure line [O IV] 
A25.9/im adopting T e ([Nelv]) and « £ ([Anv]) and the equa- 
tion of Liu et al. (2000). We chose the value from this line 
because O III lines could be affected by star light excitation 
and the abundance derived from them could be erroneous. As- 
suming the ratio of 3+ (ORLs)/0 3+ (CELs) = 10, we estimated 
the recombination contribution to [O III] A4363 less than 1 % 
of the observed value, which has a negligible effect on the 
r e ([Olll]) derivation. 

The electron temperature in BoBn 1 ranges from 9520 
([Ol]) to 14920 K ([Neiv]). Our estimated electron tem- 
peratures except for [O II] are comparable to those of Kwitter 
et al. (2003) and Kniazev et al. (2008). Their estimated tem- 
peratures are 1 2 400- 1 3 720 K for [O III] , 1 1 320- 1 1 700 K for 
[Nil], and 13 250 K for [Arm]. Note that the [Oil] electron 
temperature of 8000 K of Kwitter et al. (2003) was estimated 
adopting the [S II] density of 7100 cm" 3 . Our n e -T t plane pre- 
dicts that the [Oil] temperature is ^10000 K when adopting 
[S II] density. 

The ionic abundances derived from the CELs depend 
strongly on the electron temperature. In the case of 
[Olll]A5007, for example, only 500 K change makes a dif- 
ference of over 10% for 2+ abundance. It is therefore essen- 
tial to find the proper electron temperature for each ionized 
stage of each ion. To that end, we examined the behavior 
of the electron temperature and density as a function of LP. 
The upper panel of Fig. 8 shows that T e is increasing propor- 
tional to LP. The observed behavior of T e is consistent with the 
schematic picture of stratified physical conditions in ionized 
nebula, where the electron temperature of ions in the inner 
part should be hotter than that in the outer part. n £ is sim- 
ply monotonically increasing up to ^40 eV as I. P., except for 
[S II]. The [S II] might be emitted in high-density blobs in the 
outer nebula as we mentioned above. 

To minimize the estimated error for ionic abundances due 
to electron temperature, we have assumed a 7-zone model 
for BoBn 1 by reference to Fig. 8. Adopted T e and n £ for 
each ion are presented in Table 9. r £ ([Ol]) and n £ ([Nl]) are 
adopted for ions in zone 0, which have <10 eV. r e ([Nll]) 
and n £ ([S II]) are adopted for ions in zone 1 (LP < 1 1.3 eV). 
r £ ([N II]) and n £ ([0 II]) are for zone 2(11 .3-20 eV). T e ([S III]) 
and n £ (Cm]) are for zone 3 (20-25 eV). For zone 4 (25-41 
eV), 5 (41-63.5 eV), and zone 6 (> 63.5 eV), we adopted 
n £ ([Ar IV]) and T e ([0 III]), the averaged value from T £ ([0 III]) 
and r £ ([Ne IV]), and T £ ([Ne IV]), respectively. 

3.3.2. ORL diagnostics 



We detected a large number of optical recombination lines 
(ORL). C III.IV, 11,111,1V, N II,III, and Ne II are for the first 
time detected. To calculate ORL abundances, the electron 
temperature and density derived from ORLs are needed. We 
estimated the electron temperature using the Balmer discon- 
tinuity and He I line ratios, and the electron density using the 
Balmer decrement. The results are listed in Table 8. 

The ratio of the jump of continuum emission at the Balmer 
limit at 3646A (BJ) to a given hydrogen emission line depends 
on the electron temperature. Following Liu et al. (2001), 
we use this ratio to determine the electron temperature. This 
temperature, r £ (BJ), is used to deduce ionic abundances from 
ORLs. Defining BJ as 7(3646)-/(3681), and taking the emis- 
sivities of H I Balmer lines and H I, He I, and He II continuum 
emissivities, Liu et al. (2001) gave the following equation: 



W B,>^ 8 ( 1+ 0,59^ + 3.40^ 



/ 7(3646) -7(3681) 
v 7(H11) 



-1.5 



(5) 



[7(3646)- 7(3681 )]/7(H 11) is in units of A" 1 . T £ (BJ) is 
valid over a range from 4000 to 20 000 K. The process was 
repeated until self-consistent values for the N(He + )/N(H + ), 
N(He 2+ )/N(H + ) and T £ (BJ) were reached, we estimated T £ (BJ) 
of 8840 K. 

We estimated the He I electron temperature T^Hel) from 
the ratios of He IA7281/A6678, A7281/A5876, A6678/A4471, 
and A6678/A5876 assuming a constant electron density = 10 4 
cm" 3 , estimated from the Balmer decrement as described be- 
low. All the He I line ratios we chose here are insensitive 
to the electron density. We adopted the emissivities of He I 
from Benjamin et al. (1999). We estimated T e (He I) values as 
7340-9920 K. The r £ (Hel) from A7281/A6678 ratio appears 
to be the most reliable value because (i) He I A6678 and A728 1 
levels have the same spin as the ground state and the Case B 
recombination coefficients for these lines by Benjamin et al. 
(1999) are more reliable than the other He I A4471 and A5876; 
(ii) the effect of interstellar extinction is less due to the close 
wavelengths. We adopted 9430 K as r £ (Hel). Note that in 
Fig. 7 the electron temperatures and densities derived from 
the H I and He I are not presented. 

The intensity ratios of the high-order Balmer lines Hn (n > 
10, n: the principal quantum number of the upper level) to a 
lower Balmer line, e.g., H/3, are also sensitive to the elec- 
tron density. In Fig. 9, we plot the ratio of higher-order 
Balmer lines to H/3 with the theoretical values by Storey & 
Hummer (1995) for the cases of electron temperature of 8840 
K (= T £ (BJ)) & electron densities of 1000, 5000, 10 4 , and 
10 5 cm" 3 . This diagram indicates that the electron density 
in the ORL emitting region is between 5000 and 10 4 cm , 
which is fairly compatible with the CEL electron densities. 
Zhang et al. (2005b) estimated T £ (Hel) from the ratio of 
He I A7281/A6678 for 48 PNe and found that high-density 
blobs (10 5 — 10 6 cm" 3 ) might be present in nebula if T £ (He I) ~ 
T £ (BJ). Fig. 9 indicates that such components do not coexist 
in BoBn 1. 

3.4. Ionic abundances from CELs 

The derived CEL ionic abundances X m+ /H + are listed in Ta- 
ble 10. X" ,+ and H + are the number densities of an m times 
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FIG. 7. — Plasma diagnostic diagram. Each curve is labeled with an ID number given in Table 8. For T C ([N II]) and TJ[0 II]), we corrected for recombination 
contributions to [N II] A5755 and [O II] AA7320/30, respectively (see text). 



TABLE 8 
Plasma diagnostics. 





ID Diagnostic 


Ratio 


Result 




(1) [NelV] (A2422+A2425)/(A4715/16/25/26) 


100.78±10.84 


14 920±810 


(K) 


(2) [O III] (A4959+A5007)/(A4363) 


85.01±4.32 


13 650±290 




(3) [ArIII](A7135)/(A5192) 


85.70±35.63 


13 330±3 310 




(4) [Nem] (A15.5/im)/(A3869+A3967) 


0.20±0.()1 


13050±140 




(5) [Ne III] (A3869+A3967)/(A3344) 


333.67±14.63 


12 870±170 




(6) [Nil] (A6548+A6583)/(A5755) 


57.55±1.77 


12000±190° 




(7) [Sm](A9069)/(A6312) 


7.42±0.53 


12460±490 




(8) [Ol](A6300+A6363)/(A5577) 


68.65±10.76 


9520±550 




(9) [Oil] (A3726+A3729)/(A732()+A7330) 


10.47±0.21 


12100±180* 




(10) [S II] (A6716+A6731)/(A4069+A4076) 


13.84±5.17 


12 4 20-3590 




Average^ 




13 050 




HeI(A7281)/(A6678) 


0.21±().()1 


9430±310 




He I (A7281)/(A5876) 


0.05±0.01 


7340±110 




He I (A6678)/(A4471) 


0.83±0.02 


7400+ i070 




Hel(A6678)/(A5876) 


0.27±0.01 


9920±310 




Average 




8520 




(Balmer Jump)/(H 11) 




8840±210 


n e 


(11) [Nl](A5198)/(A5200) 


1.43±0.03 


103()±13() 


(cm' 


' 3 ) (12) [Oll](A3726)/(A3729) 


1.65±0.03 


1510±60 




(13) Cm](A1906)/(A1909) 


1.39±0.03 


3590±1000 




(14) [ArlV] (A4711)/(A4740) 


1.05±0.07 


3960±1090 




(15) [Nelll](A15.5/im)/(A36.0/mi) 


12.11±1.69 


4400+ 90io 




(16) [S II] (A6716)/(A6731) 


8.57±0.03 


5740±1310 




Average^ 




3370 




Balmer decrement 




5000-10 000 



^rom ions with I.P> 13.6 eV. 

"Corrected for recombination contribution to [NlI]A5755 (see text). 
b Corrected for recombination contribution to [O II] AA7320/30 (see text). 
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TABLE 9 

Adopted T t and n e for CEL ionic abundance calculations. 



Zone 


Ions 


r E (K) 


n e (cm 3 ) 





C ,N°,O°,Ba + 


9520 


1030 


1 


C+.S+ 


12000 


5740 


2 


+ ,N+,F + ,Fe 2+ 


12000 


1510 


3 


C 2+ ,Ne+,S 2+ ,Cl 2+ ,Xe 2+ 


12460 


3590 


4 


N 2+ ,0 2+ ,Ne 2+ Ar 2+ Ar 3+ ,S 3+ 


13650 


3960 




Cl 3+ ,F 2+ ,Fe 3+ ,Kr 3+ ,Kr 4+ 






5 


C 3+ ,N 3+ ,0 3+ 


14290+ 


3960 


6 


Ne 3+ ,Ne 4+ ,F 3+ 


14920 


3960 



The averaged value from T e ([0 III]) and 7/ e ([Ne IV]). 
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FIG. 9. — Plot of the intensity ratio of the higher order Balmer lines (Hn, 
n= 1 1-25, n: quantum number of the upper level) to H/3 (Case B assumption) 
with the theoretical intensity ratios for T € = 8840 K & and different n e 's. 



ionized ion and ionized hydrogen, respectively. To estimate 
ionic abundances, we solved level populations for a multi- 
level atomic model. In the last one of the line series of each 
ion, we present the adopted ionic abundances in bold face 
characters. These values are estimated from the line intensity- 
weighted mean or average if there are two or more available 
lines. Over 10 ionic abundances of some elements are esti- 
mated for the first time. These newly estimated ionic abun- 
dance would reduce the uncertainty of estimation of each ele- 
mental abundance, in particular, N, O, F, Ne, S, Fe, and some 
i-process elements, which are key elements to the nucleosyn- 
thesis in low-mass stars and chemical evolution in galaxies. 

Ne 2+ (zone 4 ion) and S 2+ (zone 3) abundances are derived 
from CELs seen in both the UV-optical and mid-infrared re- 
gions. CELs in the mid-infrared, namely, fine-structure lines, 
have an advantage in derivations of ionic abundances. Since 
the excitation energy of the fine-structure lines is much lower 
than that of the other transition lines, ionic abundances de- 
rived from these lines are nearly independent of the electron 
temperature or temperature fluctuation in the nebula. Note 
that these ionic abundances derived from fine-structure lines 
are almost consistent with those from other transition lines. 
This means that the adopted electron temperature and density 
for the ions in zones 3 and 4 are appropriate at least. 

Followings are short comments on derivations of ionic 
abundances. We subtract the [Olll] A2322 contamination 
from C II] A2324 using the theoretical intensity ratio [O III] 
/(A2322)/(A4363) of 0.24 and then estimate the C + abun- 
dance. The N + and + abundances are derived only from 
[N II] AA6548/83 and [O II] AA3726/29 to avoid recombina- 
tion contamination, respectively, while The Ne + abundance 
is calculated from [Nell] A12.8 /im by solving a two level 
atomic model. 

We have detected two fluorine lines [Fiv] AA3996,4060 
and estimate a F 3+ abundance of 1.50(-8). Otsuka et al. 
(2008a) detected these fluorine lines and estimated a F 3+ abun- 
dance of 5.32(-8). The F 3+ abundance discrepancy between 
Otsuka et al. (2008a) and the present work is due to dif- 
ferent adopted electron temperature and c(H/3) values. We 
have detected candidates of [Fll] AA4790,4868 and [Fill] 
AA5721/33. In the previous section, we confirmed that BoBn 
1 has no high-density components, larger than the critical 
density of these [Fill] lines. The critical density of [Fll] 
AA4790/4868 is ~2x 10 6 cm" 3 , and that of [Fin] AA5721/33 
is ^8xl0 6 cm" 3 . Therefore, the effect of collisional de- 
excitation is negligibly small. Accordingly, the ratios of 
[Fll] 7(A4790)/(A4868) and [Fill] 7(A5733)/(A5721) depend 
on their transition probabilities. When adopting the tran- 
sition probabilities by Baluja & Zeippen (1988) for [Fll] 
and Naqvi (1951) for [Fill], the theoretical intensity ratios 
of [Fll] 7(A4790)/(A4868) and [Fill] 7(A5721)/(A5733) are 
~3.2 and ~1, which are in agreement with our measurements 
(4.2±1.0 for [Fir] and 1.0±0.3 for [Fill]). Hence, these 
four emission lines can be identified as [Fll] AA4790,4868 
and [Fill] AA5721/33. The F 2+ and F 3+ abundances are es- 
timated from the each detected line by solving the statistical 
equilibrium equations for the lowest five energy levels. For 
[F III] lines the relevant collision strength has not been cal- 
culated. However, since F 2+ is isoelectronic with Ne 3+ , and 
collision strengths for the same levels along an isoelectronic 
sequence tend to vary with effective nuclear charge (Seaton 
1958). We therefore assume that the collision strengths for 
[F III] are ~22% smaller than those for [Ne IV] and estimate 
F 3+ abundances. Otsuka et al. (2008a) showed a correlation 
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between [Ne/Ar] and [F/Ar] in PNe, suggesting that Ne and F 
were synthesized in the same layer and carried to the surface 
by the third dredge-up. If this is the case, the ionic abundance 
ratios of F 2+ (LP = 35 eV) and F 3+ (62.7 eV) to F+ (17.4 eV) 
should be comparable to those of Ne 2+ (41 eV) and Ne 3+ (63.5 
eV) to Ne + (21.6 eV). Indeed, these ionic abundance ratios 
follow our prediction; F^F^F 34- -1:34:1 and Ne + :Ne 2+ :Ne 3+ 
— 1:28:1. This means that our identifications of two [Fll] and 
[Fill] lines and the estimated ionic abundances are reliable. 
So far, fluorine has found in only a handful of PNe (see Zhang 
& Liu 2005; Otsuka et al. 2008a). Among them BoBn 1 ap- 
pears to be the most F-rich PN. 

We subtract the contribution to [CI III] A8500.2 due to C III 
A8500.32 using the C 3+ ORL abundance and give upper limit 
of Cl 2+ abundance from this line. The adopted CI abun- 
dance is from [CI III] A5517 only. The Ar 3+ abundance is 
from [Ariv] AA47 17/40. We adopt the S + abundance based 
on [S II] lines except [S II] A4068, because [S II] A4068 could 
be partially blended with C III A4068. To estimate Fe 2+ and 
Fe 3+ abundances, we solved a 33 level model (from 5 Dj, to 
b 3 f 2 ) for [Fein] and a 18 level model (from 6 S 5/2 to 2 F 5/2 ) 
for [Feiv]. We adopt the transition probabilities of [Feiv] 
recommended by Froese-Fischer & Rubin (1998). For those 
not considered by Froese-Fischer & Rubin (1998), the values 
by Garstang (1958) were adopted. 

3.5. Ionic abundances of heavy elements (Z > 30) 

We have detected 10 emission line candidates of krypton 
(Kr), rubidium (Rb), xenon (Xe), and barium (Ba). Kr and Rb 
are light- s-process elements (30 ^ Z ^40, Z: atomic num- 
ber), Xe and Ba are heavy s-process (Z ^ 41). Kr has been 
detected in over 100 PNe, while the latter three ^-process ele- 
ments have been detected in only a handful of PNe (Sharpee 
et al. 2007). Selected line profiles of these candidates are pre- 
sented in Fig. 10. The Kr 3+ , Kr 4+ , Xe 2+ , and Ba + abundances 
in this object are estimated for the first time. 

We have detected two nebular lines of [Kriv] 
AA5346.7,5867.7 ( 4 S° /2 - 2 D° 5/2 and 4 S°y 2 - 2 D°y 2 , respec- 
tively). For [Kriv] A5346.7, the possibility of blending with 
C III A5345.85 (multiplet V13.01) is low, and the contribution 
to this [Kriv] line is probably negligible because other 
V13.01 Cm lines are not detected. For [Kriv] A5867.7, 
we estimated the contamination from He II A 5867.7 us- 
ing theoretical ratios of Hen /(A5867.7) to /(A5828.4), 
/(A5836.5), /(A5857.3), /(A5882.12), and /(A5896.8) given 
by Storey & Hummer (1995) assuming T e = 8840 K and 
n e = 10 4 cm" 3 , from which the contribution from Hell 
A5867.7 is estimated to be —64 %. In Fig. 10, we present 
the isolated [Kriv] A5867.7 profile. The observed ratio of 
[Kriv] /(A5867.7)//(A5346.7) (<1.8) is comparable to the 
theoretical value (—1.5) if the population at each energy 
level does not exceed the critical densities of — 1.5(+6) cm" 3 
( 2 D° 5/2 ) and -2(+7) cm" 3 ( 2 D" 3/2 ). We therefore identified 
these two lines as [Kriv] AA5346.7,5867.7. 

[Krv] A8243.4 ( 3 P 2 - l D 2 ) is likely to be blended with the 
blue wing of the Paschen series H I A8243.7 (n=3-43). Us- 
ing theoretical ratios of Hi /(A8243.7) to /(A8247.7) and 
/(A8245.6) given by Storey & Hummer (1995), we sub- 
tracted the contribution of Hi A8243.7, then estimated the 
intensity of [Krv] A8243.4. We found another nebular line 
[Krv] A6256.1 ( 3 Pi- l D 2 ). The theoretical intensity ratio of 
/(A6256.1)//(A8243.7) (—1.1) is in good agreement with ours 
(1.2). 



[Xe III] A5846.8 ( 3 P 2 - l D 2 ) appears to be blended with He II 
A5846.7. We subtracted the He II A5846.7 contribution from 
it using the theoretical ratios of /(A5846.7) to /(A5828.4), 
/(A5836.5), /(A5857.3), /(A5882.12), and /(A5896.8) given 
by Storey & Hummer (1995), then obtained an upper limit 
to the intensity of [Xe III] A5846.8. In Fig. 10, we present the 
isolated [Xe III] A5846.8 profile. 

Two Ball recombination lines AA4934,6141.7 
(6s 2 S l/2 -6p 2 P° /2 and 5d 2 D 5/2 -6p 2 P° /2 ) are detected. 
Following Sharpee et al. (2007), we estimated the Ba + 
abundances adopting transition electron temperature and 
density (zone 0). The Ba + abundances from these lines are in 
good agreement each other. 

We have detected a candidate [Rb v] A5363.6 (auroral line; 
4 Sy 2 - 2 Dy 2 ). Rb is one of the important elements as tracers 
of the neutron density. In the case of NGC 7027, Sharpee 
et al. (2007) argued the possibility that this line is Oil 
A5363.8 (4fF 2 [4]° /2 -4d' 2 F 7/2 ). They also suggested that 
the intensity of Oil A5363.8 is comparable to Oil A4609.4 
(3d 2 D 5/2 -4fF 2 [4] y2 ), arising from the lower level of Oil 

A5363.8. In BoBn 1, based on the ORL 2+ abundance of 
1.45(-4) from the 3d-4f Oil lines (see next Section) the ex- 
pected intensity of Oil A4609.4 is -6.8(-5)/(H/3), which is 
lower than the observed intensity of the [Rb v] A5363.6 can- 
didate. No 4fF 2 [4]°-4d' 2 F Oil lines are detected in BoBn 
1. Therefore, we consider that the detected line is [Rbv] 
A5363.6. Since there are no available collision strengths for 
this line at present, we do not estimate a Rb 4+ abundance. 

3.6. Ionic abundances from ORLs 

We have detected many optical recombination lines (ORLs) 
of helium, carbon, nitrogen, oxygen, and neon. To our knowl- 
edge, the nitrogen, oxygen, and neon ORLs are detected for 
the first time from this PN. These lines provide us with a 
new independent method to derive chemical abundances for 
BoBn 1. The recombination coefficient depends weakly on 

-1/2 

the electron temperature (oc T e ' ). The ionic abundances are, 
therefore, insignificantly affected by small-scale fluctuations 
of electron temperature. This is the most important advantage 
of this determination method. The ionic abundances from re- 
combination lines are robust against uncertainty in electron 
temperature estimation. 
The ORL ionic abundances X m+ /H + are derived from 

X m+ _ a(H/3) A(X m+ ) 7(X m+ ) 

IF" ~ a(X m+ ) A(H/3) 7(H/3) ' () 

where a(X m+ ) is the recombination coefficient for the ion 
X m+ . For calculating ORL ionic abundances, we adopted T t of 
8800 K and « £ of 10 4 cm" 3 from the hydrogen recombination 
spectrum. 

Effective recombination coefficients for the lines' parent 
multiplets were taken from the references listed in Table 1 1 . 
The recombination coefficients for each multiplet at a given 
electron density were calculated by fitting the polynomial 
functions of T t . The recombination coefficient of each line 
was obtained by a branching ratio, B(A,), which is the ratio of 
the recombination coefficient of the target line, a(A,) to the 
total recombination coefficient, 5Z,-a(Ai) in a multiplet line. 
To calculate the branching ratio, we referred to Wiese et al. 
(1996) except for O II 3p-3d and 3d-4f transitions and Ne II. 
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FIG. 10.— Detected fluorine lines and candidates of ^-process [KrIV]A5867.7, [KrV]A6256.1, [Rb V]A5363.6, [Xe III]A5846.7, and BaIlA6141.7. The thin 
lines indicate the fitted Gaussian profiles to each emission line. The isolated [Kr IV] A5867.7 and [Xe III] A5846.7 line-profiles after smoothing are also presented 
(see text for detail). 



For O II 3p-3d and 3d-4f transition lines, the branching ratios 
were provided by Liu et al. (1995) based on intermediate cou- 
pling. For Ne II, Kisielius et al. (1998) provided the branching 
ratios based on LS-coupling. 

The estimated ORL ionic abundances are listed in Tables 
12 and 13. In general, a Case B assumption applies to the 
lines from levels having the same spin as the ground state, and 
a Case A assumption applies to lines of other multiplicities. 
In the last one of the line series of each ion, we present the 
adopted ionic abundance and the error, which are estimated 
from the line intensity-weighted mean. 

3.6.1. Helium 

The He + abundances are estimated using electron density 
insensitive five He I lines to reduce intensity enhancement by 
collisional excitation from the He 2s 3 S level. The collisional 
excitation from the He 2s 3 S level enhances mainly the in- 
tensity of the triplet He I lines. We removed this contributions 
(1.4% for He I A4387; up to 7.4% for He I A5876) from the 
observed line intensities using the formulae given by King- 
don & Ferland (1995). 

The He 2+ abundance is estimated from He II A4686. Kni- 
azev et al. (2008) estimated He + = 8.52(-2) and He 2+ = 1.53(- 
2), which are close to or slightly smaller than our values. 

3.6.2. Carbon 

We observed C II lines which arose from different transi- 
tions. The ground state of C II line is a doublet (2p 2 Pq). The 
3d-4f (multiplet V6), 4d-6f (V16.04), 4f-6g (V17.04), and 4f- 
7g (V17.06) lines, which have higher angular momentum as 
upper levels, are unaffected by both resonance fluorescence 
by starlight and recombination from excited 2 S and 2 D terms. 
Among these high angular momentum lines, the V6 lines are 
the most case-insensitive and reliable. Comparison of the C 2+ 
abundance derived from C II A4267 with that of the other C II 
lines indicates that the observed C II lines are not populated 
by the intensity enhancement mechanisms discussed above. 
Therefore we can safely use all the C II lines for the estima- 
tion of C 2+ abundance. 



All the observed C III lines are triplets. Since the ground 
state of C III is singlet (2s 2 1 S), we adopted Case A assump- 
tion. Unlike the case of C II, C III lines are relatively case in- 
sensitive. Our estimated C 2+ and C 3+ abundances (Table 12) 
are in good agreement with Kanizev et al. (2008); their C 2+ 
and C 3+ are 7.78(^0 and 5. 62(^4), respectively. 

We estimate the C 4+ abundance using multiplet V8 and 
V8.01 lines. Interestingly, we observed Civ A5811. Civ 
A5801/1 1 has been detected in PNe with Wolf Rayet type cen- 
tral stars, suggesting that the central star is very active. In the 
case of BoBn 1 C IV lines might be be nebular origin rather 
than the central star origin, because the 2V exp of C IV A58 1 1 is 
14.3 km s" 1 comparable with the value in close LP. ions such 
as [Ne IV] and [Fiv] (see Table 6). 

3.6.3. Nitrogen 

All of the observed N II lines are triplets. Since the ground 
level of N II is a triplet (2p 2 3 P), we adopted Case B assump- 
tion. The Nil resonance line 2p 2 3 P-2p4s 3 Pi A508.668 A 
can be enhanced by the He I resonance line 1 s 2 1 S-ls&p 1 Pq 
A508.697 A. The cascade transition from 2p4s 3 Pi can en- 
hance the line intensity of the multiplet V3 lines. But, this 
transition cannot enhance the line intensity of 3f-4d transition 
(multiplet V43b, V48a, V50a, and V55a) due to the lack of 
a direct resonance or cascade excitation path. Comparison of 
N 2+ abundances derived from the 3f-4d with those from the 
V3 lines implies that the fluorescence is negligible in BoBn 1 . 

The multiplet VI, V2, and V17 N III lines are observed. We 
adopted Case B assumption except for the V17 multiplet. For 
the V17 line, we adopted Case A assumption. The intensity 
of the resonance N III line A374.36 A (2p 2 P°-3d 2 D) may 
be enhanced by O m resonance at 374.1 1 A (2p 2 3 P-3s 3 P°). 
The line intensity of the multiplet VI and V2 lines might be 
enhanced by the Olll lines. The multiplet V17 line (4f-5g 
transition) does not appear to be enhanced. Therefore, we 
adopt the N 3+ abundance from this line. 

3.6.4. Oxygen 
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TABLE 10 
Ionic abundances from CELs. 



TABLE 10 
Continued. 



X m+ A lab 


/(Alab) 


T e n t 


X m+ /H+ 


X" 1+ A, ab 


/(Alab) 


T e n e 


X m+ /H+ 


(A//im) 


[7(H/3)=100] 


(K) (cm" 3 ) 




(A//im) 


[7(H/3)=100] 


(K) (cm" 3 ) 





c + 

c 2+ 



c 3+ 



N+ 



N 2+ 
N 3+ 

o° 



+ 



2+ 



o 3 + 

F+ 



p2+ 



8727.12 7.93(-2) ± 4.73(-3) 
2324 3.51(+1)± 4.46(0) 
1906 8.39(+2) ± 1.27(+1) 
1908 6.03(+2) ± 1.03(+1) 



1548 
1551 



1.05(+3)±2.02(+l) 
5.19(+2) ± 1.50(+1) 



9520 1030 
12000 5740 
12460 3590 



14290 3960 



N° 5197.90 2.73(-l) ± 7.90(-3) 9520 1030 
5200.26 1.91(-1) ± 5.47(-3) 



5754.64 
6548.04 
6583.46 

1750 
1485 
5577.34 
6300.30 
6363.78 

3726.03 
3728.81 
7319 
7330 

4363.21 
4931.23 
4958.91 
5006.84 

25.9 
4789.45 
4868.99 



1.23(0) ± 1.40(-2) 
1.56(+1) ± 8.34(-l) 
5.04(+l) ± 1.28(0) 

4.81(+1) ± 1.32(+1) 
4.58(+l) ± 1.78(+1) 
1.70(-2) ±2.64(-3) 
8.72(-l) ± 1.55(-2) 
2.91(-1) ± 9.91(-3) 

1.09(+1) ± 6.85(-2) 
6.61(0) ± 9.99(-2) 
1.02(0) ± 2.00(-2) 
7.81(-1) ± 1.38(-2) 

5.57(0) ± 9.73(-2) 
4.36(-2) ± 3.88(-3) 

1.22(+2)± 5.96(0) 
3.51(+2)±2.18(+1) 

1.25(+1) ± 1.36(-1) 
5.58(-2) ± 5.20(-3) 
1.34(-2)±3.00(-3) 



12000 1510 



13650 3960 
14290 3960 
9520 1030 



12000 1510 



13650 3960 



14290 3960 
12000 1510 



5721.20 2.70(-2) ± 2.93(-3) 13650 3960 
5733.05 2.68(-2) ± 5.43(-3) 

3996.92 4.09(-2) ± 2.53(-3) 14920 3960 
4059.90 1.19(-l)±3.72(-3) 



Ne+ 
Ne 2 " 



12460 3590 
13650 3960 



12.8 2.49(0) ± 8. 17(-2) 

3342.42 8.47(-l) ± 2.75(-2) 

3868.77 2.17(+2) ± 1.05(+1) 

3967.46 6.39(+l) ±4.13(-1) 

4011.60 1.48(-2)±3.92(-3) 

15.6 1.61(+2) ± 1.30(0) 

36 1.33(+1)± 1.84(0) 



Ne 3+ 2423.50 1.71(+1) ± 1.78(0) 14920 3960 



4714 


25 


5.52(-2) 


± 


4715 


80 


1.87(-2) 


± 


4724 


15 


5.08(-2) 


± 


4725 


62 


4.52(-2) 


± 


Ne 4+ 3345 


83 


3.22(-l) 


± 


3425 


87 


8.71(-1) 


± 


S+ 4068 


60 


3.95(-l) 


± 


4076 


35 


2.45(-2) 


± 


6716 


44 


1.23(-1) 


± 


6730 


81 


2.16(-1) 


± 



-3) 



14920 3960 



12000 5740 



4.74(-7) 
2.30(-5) 

7.71(-4) 
7.70(-4) 
7.71M) 
2.42(-4) 
2.36M) 
2.40(-4) 
4.81(-7) 
4.82(-7) 
4.82(-7) 
7.20(-6) 
5.84(-6) 
6.40(-6) 
6.27(-6) 
6.24(-5) 
3.83(-5) 
2.06(-6) 
2.04(-6) 
2.13(-6) 
2.06(-6) 
4.00(-6) 
4.03(-6) 
7.13(-6) 
6.74(-6) 
4.01(-6) 
4.76(-5) 
4.37(-5) 
4.80(-5) 
4.76(-5) 
4.77(-S) 
3.41(-6) 
2.16(-8) 
1.66(-8) 
1.98(-8) 
6.59(-7) 
6.70(-7) 
6.65(-7) 
1.47(-8) 
1.51 (-8) 
1.50(-8) 
2.97(-6) 
6.40(-5) 
8.41(-5) 
5.94{-5) 
9.71 (-5) 
9.15(-5) 
9.07(-5) 
8.32(-5) 
3.97(-6) 
4.35(-6) 
5.05(-6) 
3.60(-6) 
3.43(-6) 
3.97(-6) 
1.99(-7) 
1.97(-7) 
1.98(-7) 
4.05(-8) 
7.44(-9) 
1.03(-8) 
1.03(-8) 
1.03(-8) 



± 9.34(-8 
± 3.48(-6 

± 1.59(--4 
± 1.59(-4; 
± 1.59(-4 
± 4.47(-5 
± 4.38(-5 
± 4.44(-5; 
± 1.01(-7 
± 9.67(-8; 
± 9.90(-8 
± 4.74(-7 
± 3.75(-7 
± 2.79(-7 
± 3.02(-7; 
± 1.88(-5 
± 1.65(-5 
± 7.06(-7 
± 3.97(-7 
± 4.20(-7 
± 4.03(-7; 
± 2.23(-7 
± 2.42(-7 
± 5.67(-7 
± 5.30(-7 
± 2.30(-7 
± 4.87(-& 
± 4.58(-& 
± 3.52(-& 
± 3.94(-6; 
± 3.83(-6 
± 8.21(-8 
± 2.23(-9 
± 3.79(-9 
± 2.74(-9 
± 1.03(-7 
± 1.55(-7 
± 1.29(-7; 
± 2.21(-9 
±2.12(-9 
± 2.14(-9 

± i.i4(-7; 

± 8.10(-6; 
± 6.75(-6; 
± 3.82(-6; 
± 2.65(-5 
± 1.27(-6; 
± 1.26(-5 
± 4.33(-6 
± 8.99(-7; 
± 1.23(-7 
± 1.52(-6; 
± 1.01(-6 
± 9.78(-7 
± 9.01(-7; 
± 3.40(-8 
±3.15(-8 
± 3.22(-8 
± 1.96(-9 
± 2.79(-9 
±5.12(-10 
± 4.01(-10; 
± 4.41(-10 



6312.10 4.77(-2) ± 4.75(-3) 12460 3590 



9068.60 
18.7 



3.78(-l) ± 9.97(-3) 
6.92(-l) ± 4.67(-2) 



S 3+ 
Cl 2+ 

Cl 3+ 
Ar 2+ 



Ar 3+ 



Fe 2+ 

Fe 3+ 
Kr 3+ 

Kr 4 ' 



10.5 1.92(0) ± 5.03(-2) 
5517.66 1.81(-2) ± 2.82(-3) 
8500.20 <1.81(-3) 
8046.30 2.05(-2) ± 2.60(-3) 
5191.82 3.90(-3) ± 1.61(-3) 
7135.80 2.73(-l) ± 1.14(-2) 
7751.10 6.11(-2)±2.51(-3) 



13650 3960 

12460 3590 

13650 3960 

13650 3960 



4711 


37 


9.40(-2) 


± 


4740 


17 


8.99(-2) 


± 


7170 


50 


6.53(-3) 


± 


7262 


70 


4.58(-3) 


± 


4881 


00 


2.14(-2) 


± 


5270 


40 


2.19(-2) 


± 



13650 3960 



12000 1510 



6740.63 1.58(-2) ± 4.81(-3) 
5346.02 4.56(-3) ± 2.02(-3) 
5867.70 <8.20(-3) 
6256.06 <6.27(-3)) 
8243.39 <5.12(-3) 



Xe 2+ 
Ba + 



5846.66 <1.56(-3) 
4934.08 6.42(-3)± 1.80(-3) 
6141.70 3.43(-3) ± 7.66(^1) 



13650 3960 
13650 3960 

13650 3960 



12460 3590 
9550 1030 



6.87(-8) ± 1.07 (-8) 
7.34(-8) ± 4.79(-9) 
6.81(-8) ± 4.82(-9) 
6.99(-8) ± S.06(-9) 
5.28(-8) ± l.Sl(-9) 
1.36(-9) ± 2.42(-10) 

<2.79(-9) 
7.82(-10) ± 1.04(-10) 
1.23(-8) ± 5.17(-9) 
1.30(-8) ± 7.41(-10) 
1.22(-8) ± 6.83(-10) 
1.29(-8) ± 7.30(-10) 
7.53(-9) ± 5.46(-10) 
7.58(-9) ± 4.60(-10) 
4.32(-8) ±6.13(-9) 
3.52(-8) ± 3.08(-8) 
7.56(-9) ± 5.04(-10) 
1.16(-8)±2.80(-9) 
1.04(-8) ± 1.79(-9) 
1.10(-8) ± 2.29(-9) 
1.02(-7) ± 3.26(-8) 
1.41(-10) ± 6.31(-11) 
<1.89(-10) 
<4.09(-10) 
<3.45(-10) 
<3.77(-10) 
<2.30(-ll) 
1.90(-10)±6.45(-11) 
2.12(-10) ±6.42(-ll) 
1.98(-10) ± 6.44(-ll) 



TABLE 11 

Effective recombination coefficient References. 



Line 


Transition 


References 


Hi 


All 


(1),(2) 


He I 


Singlet 


(3) 




Triplet 


(4) 


Hell 


3-4, 4-6 


(4) 


CII 


3d-4f, 3s-3p, 3p-3d 


(5) 




4f-7g, 4d-6f, 4f-6g 




Cm 


3s-3p 


(6) 




4f-5g 


(4) 


Civ 


2p-2s, 5fg-6gh 


(4) 


Nil 


3s-3p, 3p-3d 


(7) 




3d-4f 


(8) 


Nin- 


3s-3p, 3p-3d, 4f-5g 


(4) 


On 


3s-3p 


(9) 




3p-3d, 3d-4f 


(10) 


Om 


3s-3p 


(4) 


Nell 


3s-3p, 3p -3d , 3s -3p 


(11) 



REFERENCES. — (1) Aller (1984). (2) Storey & Hummer (1995). (3) 
Benjamin et al. (1999). (4) Pequignot et al. (1991). (5) Davey et al. (2000). 
(6) Nussbaumer & Storey (1984). (7) Kisielius & Storey (2002). (8) Escalate 
& Victor (1990). (9) Storey (1994). (10) Liu et al. (1995). (11) Kisielius et 
al. (1998). 

We observed O II doublet (3d-4f) and quadruplet lines (mul- 
tiplet VI, V4, V10, VI 9). Most of the V 1 lines and all of the 
V 2 lines are observed. Since the ground level of O II is a 
quadruplet, we adopted Case A for the doublet lines and Case 
B for quadruplet lines. It seems that the multiplet VI and V10 
lines give the most reliable value. 



TABLE 13 

N, O, AND NE IONIC ABUNDANCES FROM ORLS. 



Multi. Ai ab /(Alab) 

(A) [7(H/3)=100] N 2+ /H+ 



16 Otsuka et al. 



TABLE 12 

HE AND C IONIC ABUNDANCES FROM ORLS. 



Multi. 


-Mab 


/(Alab) 






(A) 


[7(H/3)=100] 


He+/H+ 


Vll 


5876.62 


18.1 ±0.14 


9.93(-2) ± 1.26(-3) 


V14 


4471.47 


4.82 ± 0.04 


9.63(-2) ± 2.89(-3) 


V46 


6678.15 


4.02 ±0.11 


9.59(-2) ± 3.80(-3) 


V48 


4921.93 


1.32 ±0.02 


9.66(-2) ± 3.02(-3) 


V51 


4387.93 


0.59 ± 0.01 


9.54(-2) ± 3.57(-3) 




Adopted 




9.81(-2) ± 2.01(-3) 



He 2+ /H+ 



3.4 4685.68 24.8 ± 0.79 2.03(-2) ± 6.47M) 

C 2+ /H + 

V2 6578.05 3.75(-l) ± 6.87(-3) 7.30M) ± 3.50(-5) 
V6 4267.15 7.90(-l) ± 4.42(-2) 7.55M) ± 5.02(-5) 
V16.04 6151.27 3.81(-2) ± 3.40(-3) 8.74(-4) ± 8.26(^5) 
V17.04 6461.95 7.80(-2) ± 6.77(-3) 1.2A(-A) ± 6.97(-5) 
V17.06 5342.43 5.39(-2) ± 3.94(-3) ± 8.16(-5) 

Adopted 7.58(-4) ± 4.92(-5) 

C 3+ /H + 

VI 4647.42 4.41(-1) ± 3.82(-3) 7.60M) ± 2.15(-5) 
VI 4650.25 2.61(-1) ± 4.77(-3) 1A9(-A) ± 2.44(-5) 
V16 4067.87 2.49(-l) ± 8.01(-3) 6.11(-4) ± 2.80(-5) 
V16 4070.20 4.07(-l) ± 1.10(-2) 5.5A(-A) ± 2.36(-5) 
V18 4186.90 3.46(-l)± 5.45(-3) 5.80M)±2.11(-5) 
V43 8196.50 4.39(-l) ± 8.34(-3) 5.66M) ± 2.17(-5) 

Adopted 5.74(-4) ± 2.32(-5) 

C 4+ /H + 

V8 4658.64 1.19(-1) ± 1.91(-2) 2.69(-5) ± 4.41(-6) 
V8.01 7725.90 3.11(-2) ± 1.52(-3) 1.49(-5) ± 8.81(-7) 
Adopted 2.69(-5) ± 4.41(-6) 

A number of O III lines are observed. We consider Case B 
for the triplet lines (multiplet V2) and Case A for the singlet 
line (multiplet V5). There is a possibility that the multiplet V2 
lines would be excited by the Bowen fluorescence mechanism 
or by the charge exchange of 3+ and H° instead of recombi- 
nation and the multiplet V5 line could be excited by charge 
exchange. Therefore, we did not use 3+ abundances in the 
estimation of a total oxygen abundance from ORLs. 

3.6.5. Neon 

The observed Ne II lines are doublet (multiplets V9 and 
V21) and quartet lines (VI and V2). We considered Case B 
for the doublet lines and Case A for the quartet lines. The mul- 
tiplet VI and V2 lines are insensitive to the case assumption 
and are pure recombination lines (Grandi 1976). Therefore, 
we adopted the Ne 2+ abundance derived from multiplet V 1 
and V2 lines. 

3.7. Ionization Correction 

If the ionic abundances in all ionization stages are known, 
an elemental abundance will be simply the sum of its ionic 
abundances. Actually, it is, however, impossible to probe 
all of the ionization stages of an element using UV to mid- 
infrared spectra. To estimate elemental abundances, we must 
correct for unobserved ionic abundances. This correction 
was performed using ionization correction factors, ICF(X). 
ICFs(X) for each element are listed in Table 14. 

3.7.1. Helium, Carbon, Nitrogen, Oxygen and Neon 

The He abundance is the sum of He + and He 2+ . 

The C abundance is the sum of C + , C 2+ , C 3+ , and C 4+ . For 
the C abundance derived from ORLs, we corrected for un- 
seen C + assuming (C + /C)orl s = (N + /N) C el s - For the C abun- 



V3 5710.76 5.15(-3)±4.49(-3) 1.26(^1) ± 1.10(-4; 
V3 5685.26 2.52(-2) ± 2.36(-3) 6.62(-4) ± 6.48(-5; 
V3 5679.56 1.62(-2) ± 4.31(-3) 6.76(-5) ± 1.81(-5] 
V19 5001.47° 1.79(-2)±3.19(-3) 4.75(-5) ± 8.61(-6; 
V43b 4171.61 1.15(-2)±2.07(-3) 1.63(-4) ± 2.99(-5; 
V48a 4247.22 2.42(-2) ± 5.05(-3) 1.14(^1) ± 2.41(-5; 
V50a 4179.67 1.36(-2) ± 7.54(-3) 3.43(-4) ± 1.91(-4; 
V55a 4442.02 1.27(-2) ± 3.88(-3) 3.62(-^) ± l.ll(-4; 

Adopted 2.62M) ± 5.99(-5; 

N 3+ /H+ 

VI 4097.35 5.00(-l)±2.47(-2) 1.39(-3) ± 7.97(-5; 
VI 4103.39 3.14(-l)±3.72(-2) 1.75(-3) ± 2.13(^4; 
V2 4634.12 1.67(-l)±5.69(-3) 1.31(-4) ± 5.98(-6; 
V2 4640.64 3.22(-l) ± 3.62(-3) 1.41 (-4) ± 4.55(-6; 
V2 4641.85 4.88(-2) ± 5.97(-3) 1.92(^1) ± 2.42(-5; 
V17 4379.11 5.97(-2) ± 5.13(-3) 2.56(-5) ± 2.36(-6; 

Adopted 2.56(-5) ± 2.36(-6; 

Q 2+ /H+ 

VI 4638.86 1.40(-2) ± 1.21(-3) 1.27(^1) ± 9.78(-6; 
VI 4641.81 3.46(-2) ± 9.54(-3) 1.30(^1) ± 3.75(-5; 
VI 4649.13 1.86(-2) ± 1.43(-3) 3.95(-5) ± 2.38(-6; 
VI 4650.84 2.82(-2) ± 2.48(-3) 2.72(-4) ± 2.10(-5; 
VI 4661.63 2.32(-2) ± 4.71(-3) 1.85(^1) ± 3. 81(-5; 
VI 4673.73 1.85(-2)±9.08(-3) 9.87(-4) ± 4.68(-4; 
VI 4676.23 7.65(-3)±4.19(-3) 8.01(-5) ± 4.40(-5; 
V4 6721.39 2.99(-3) ± 3.24(-4) 5.13(^1) ± 5.73(-5; 
V10 4069.62 1.03(-2)±3.01(-3) 1.06(-4) ± 3.12(-5; 
V10 4069.88 1.60(-2) ± 1.60(-2) 1.03(-4) ± 3.04(-5; 
V19 4153.30 1.68(-2)±3.11(-3) 2.19(-A) ± 4.11(-5; 
3d-4f 4089.29 1.43(-2) ± 5.64(-3) 1.30(^1) ± 5.1 4(-5; 
3d-4f 4292.21* 1.76(-2) ± 5.00(-3) 6.28(^) ± 1.79(^1; 

Adopted 1.45M) ± 2.32(-5; 

Q 3+ /H + 

V2 3754.70 1.71( 1) ± 6.17(-3) 3.31H1) ± 1.54(-5; 
V2 3757.21 8.29(-2) ± 7.35(-3) 3.62(^1) ± 3.38(-5; 
V2 3759.88 6.09(-l) ± 1.97(-2) 6.47(-4) ± 2.82(-5; 
V2 3791.27 6.71(-2) ± 5.85(-3) 4.29(-4) ± 3.95(-5; 
V5 5592.37 1.07(-2) ± 1.96(-3) 4.42(-4) ± 8.19(-5; 

Ne 2+ /H+ 

VI 3694.21 4.14(-2) ± 8.49(-3) 1.28(^1) ± 2.71(-5; 
V2 3334.87 1.04(-1) ± 1.93(-2) 1.60(^1) ± 3.08(-5; 
V9 3568.50 6.95(-2) ± 7.83(-3) 2.24(-3) ± 2.62(^T 
V21 3453.07 1.25(-2) ± 3.82(-3) 4.94(-^) ± 1.53(-4; 

Adopted 1.51M) ± 2.98(-5; 

"blending line (A 5001.12,5001.47 lines). 

b blending line (A 4291.26,4291.86,4292.21,4292.98). 



dance from CELs, we corrected for C 4+ assuming (C 4+ /C)cel s 

= (C 4+ /C)oRLs- 

The N abundance is the sum of N + , N 2+ , and N 3+ . For the 
ORL N abundance, we corrected for N + assuming (N + /N)orl s 

= (N + /N)c E Ls. 

The O abundance is the sum of 0\ 2+ , and 3+ . For the 
ORL O abundance, we used only 2+ because most of the 
O III lines are not pure recombination lines. We assumed 

(0 2+ /0) RLs = (0 2+ /0)cELs. 

The Ne abundance is the sum of Ne + , Ne 2+ , Ne 3+ , and Ne 4+ . 
For the ORL Ne abundance, we corrected for the unseen Ne + , 
Ne 3+ and Ne 4+ assuming (Ne 2+ /Ne)o RLs = (Ne 2+ /Ne) CE Ls- 

3.7.2. Other elements 
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TABLE 14 

Adopted ionization correction factors (ICFs). 



4 - 



X Line ICF(X) 



He ORLs 1 
C CELs 1 



( N V 

V N 2 ++N J + J 

1 

orls 



ORLs 
CELs 1 



O CELs 1 
ORLs 



,o°+) f 



X/H 

He++He 2+ 

c++c 2+ +c 3+ 

ICF(C)(C 2+ + C 3+ + C 4+ ) 

N++N 2+ +N 3+ 
ICF(N)(N 2+ +N 3+ ) 

o++o 2+ +o 3+ 

ICF(0)0 2+ 



CELs ( wZ&Z& y ICF(F)(F+ + F 2 + + F 3 



Ne CELs 1 
ORLs 

S CELs 

CI CELs 

Ar CELs 

Fe CELs 

Kr CELs 

Xe CELs 



_Ne 

. Ne : 



I N \t 
V N++N 2 + / 



Ne++Ne 2+ +Ne 3+ +Ne 4+ 
ICF(Ne)Ne 2+ 

ICF(S)(S+ + S 2+ + S 3+ ) 



ICF(C1)(C1 2 + + C1 3 +) 
(n^nw)' ICF(Ar)(Ar 2 + + Ar 3 +) 



Ne 

, Ne+tNe^+Ne 4 * , 

o y 

o++o 2 + / 



Cl 3 + 



s 



Ba CELs* 1 



ICF(Fe)(Fe 2+ + Fe 3+ ) 

ICF(Kr)Kr 3+ +Kr 4+ 

ICF(Xe)Xe 2+ 

Ba + 



Ionic and elemental abundances derived from CELs. 



The value is the lower limit (see text). 



Assuming that the F abundance is the sum of F 1 ", F 2+ , F 3+ , 
and F* + , we corrected for unseen F 4+ using the CEL Ne abun- 
dance. The S abundance is the sum of S + , S 2+ , and S 3+ , 
and S 4+ . Unseen S 4+ was corrected for assuming S 4+ /S = 
(N 3+ /N)cels- We assume that the CI abundance is the sum 
of Cl + , Cl 2+ , Cl 3+ , and Cl 4+ . The unseen Cl + and Cl 4+ are 
corrected for assuming C1/(C1 + +C1 4+ ) = 0/(0 + +0 3+ ) CE Ls- For 
Ar, its abundance is assumed to be the sum of Ar 2+ , Ar 3+ , and 
Ar 4+ , and unseen Ar 4+ was corrected for assuming (Ar 4+ /Ar) 
= (Ne 4+ /Ne)cELs- For Fe, we assume that its abundance is the 
sum of Fe 2+ , Fe 3+ , and Fe 4+ . The unseen Fe 4+ was corrected 
for assuming (Fe 4+ /Fe) = (0 3+ /0) C el s - 

We assume that the Kr abundance is the sum of Kr 2+ , 
Kr 3+ , and Kr 4+ , the unseen Kr was corrected for assuming 
Kr 2+ /Kr 3+ = C1 2+ /C1 3+ . We assume that the Xe abundance is 
the sum of Xe + , Xe 2+ , Xe 3+ , and Xe 4+ . The Xe ionic abun- 
dances except Xe 2+ were corrected for assuming Xe 2+ /Xe = 
S 2+ /S. We give the lower limit of the Ba abundance, which is 
equal to the Ba + abundance (IP = 5.2 eV), since we could not 
detect higher excited Ba lines with IP > 13.5 eV. It should 
take care in handling the Ba abundance. 

3.8. Elemental abundances 

The resultant elemental abundances are presented in Ta- 
ble 15. We recognized that BoBnl is a C-, N-, and Ne-rich 
PN: the [C/O], [N/O], and [Ne/O] abundances from ORLs are 
+ 1.23, +1.12, +0.81. The ratios derived from CELs are +1.58, 
+ 1.10, and +1.04, respectively. Comparing the C, N, O, Ne 
ORL and CEL abundances, ORLs might be emitted from 0-, 
Ne-rich region. The ORL C, N, O, Ne abundances are larger 



&_ 

Q 

< 2 




10" 1.1x10" 1.2x10" 1.3x10" 1.4x10" 1.5x10" 
CELs Electron Temperature (T c ) [K] 

FIG. 1 1 . — The abundance discrepancy factor (ADF) vs. the electron tem- 
perature from CELs. The filled circles are estimated values when adopting 
the T 6 and n e values listed in Table 9. 

by 0.14-0.49 dex than the CEL abundances. We need to look 
for reasons for the abundance discrepancy. 

4. DISCUSSION 

First, in this section, we will discuss the abundance discrep- 
ancies between CELs and ORLs using three models (Section 
4.1). Second, we will compare elemental abundances esti- 
mated by us with others (Section 4.2). Third, we will build a 
photo-ionization model to derive the parameters of the central 
star, ionized nebular gas, and dust (Section 4.3). Next, the em- 
pirically derived elemental abundances will be compared with 
theoretical nucleosynthesis model predictions for low- to in- 
termediate mass stars (Section 4.4). Finally, we will guess 
the evolutionary status or provide a presumable evolutionary 
scenario for BoBn 1 (Section 4.5). 

4.1. The abundance discrepancy between CELs and ORLs 

We derived ionic and elemental abundances using CELs 
and ORLs and found somewhat large abundance discrepan- 
cies between them. So far, abundance discrepancies have 
been found in about 90 Galactic disk PNe, 3 Magellanic PNe 
(Tsamis et al. 2003, 2004; Liu et al. 2004; Robertson-Tessi 
& Garnett 2005; Wesson et al. 2005, etc.), and 1 Halo PN 
(DdDm 1; Otsuka et al. 2009). We define the ionic abundance 
discrepancy factor ADF as the ratio of the ORL to the UV or 
optical CEL abundances. In BoBn 1, the ADFs are 0.98±0.21 
for C 2+ , 2.39±0.45 for C 3+ , 4.2 1 ± 1 .59 for N 2+ , 0.67±0.29 for 
N 3+ , 3. 05 ±0.54 for 2+ , and 1 .82±0.39 for Ne 2+ , respectively 

Up to now, three models have been proposed to explain 
abundance discrepancies in PNe: temperature fluctuations, 
high density components, and hydrogen-deficient cold com- 
ponents. We examine what can cause the abundance discrep- 
ancies in BoBn 1 using these models. 

4.1.1. Temperature fluctuations 

The emissivities of the CELs increase exponentially as the 
electron temperature becomes higher. The electron tempera- 
ture derived from the CELs, r c (CELs) will be indicative of the 
hot region nearby the radiation source. If we adopt r e (CELs) 
for abundance estimations using the CELs, the ionic abun- 
dances might be underestimated. 

Peimbert (1967) considered the effect of electron tempera- 
ture fluctuation in a nebula, which sometimes gave high elec- 
tron temperature, on the determinations of the ionic abun- 
dances derived from CELs. For example, Torres-Peimbert 
et al. (1980) characterized the electron temperature fluctu- 
ations in term of t 1 as the cause of the abundance discrep- 
ancy between CELs and ORLs. Assuming the validity of the 
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TABLE 15 

The elemental abundances derived from CELs and ORLs in the case of no temperature fluctuation. 



CELs ORLs CELs ORLs CELs ORLs 



X/H 



log(X/H)+12 a 



[X/H] 6 



He 
C 
N 
O 
F 
Ne 
S 
CI 
Ar 
Fe 
Kr 
Xe 



1.05(-3) ± 1.93(-4) 
1.07 HI) ± 2.50(-5) 
5.51(-5) ± 3.84(-6) 
7.01(-7) ± 1.38(-7) 
9.04(-5) ± 4.42(-6) 
2.07(-7) ± 7.53(-8) 
2.47(-9) ± 4.02(-10) 
2.13(-8) ± 1.75(-9) 
1.21(-7) ± 3.69(-8) 
<7.63(-10) 
<9.33(-ll) 



1.18(-1) ± 2.12(-3) 
1.44(-3) ± 4.96(-4) 
3.06(-4) ± 1.22(-5) 
1.68(-4) ± 3.22(-5) 

1.64(-4) ± 3.44(-5) 



Ba 1.98(-10) ±6.44(-ll) 



9.02±0.08 
8.03±0.10 
7.74±0.03 
5.85±0.09 
7.96±0.02 
5.32±0.17 
3.39±0.07 
4.33±0.04 
5.08±0.14 
<2.88 
<1.97 
2.30±0.15 



11.07±0.01 
9.16±0.16 
8.49±0.18 
8.23±0.08 

8.22±0.09 



+0 

+0.63±0.09 +0 
+0.1 5 ±0.15 +0 
-0.95±0.06 -0 
+1.39±0.11 
+0.09±0.10 +0 
-1.87±0.17 
-1.94±0.09 
-2.22±0.09 
-2.39±0.14 

<-0.48 

<-0.27 
+0.12±0.15 



"The number density of the hydrogen is 12. 
b Solar abundances are taken from Lodders (2003). 

temperature fluctuation paradigm, the comparison of the ionic 
abundances derived from CELs and ORLs may provide an es- 
timation of f 2 . 

The relation between ADFs and r e (CELs) is presented in 
Fig. 11. We recognize that ADFs would approach to 1 
if r £ (CELs) for each zone dropped by >1000 K. Using the 
formulations for temperature fluctuations given by Peimbert 
(1967), we have estimated the f 2 parameter and the mean elec- 
tron temperatures Tq for each zone . The resultant f 2 and Tq 
are listed in Table 16. The derived f 2 = 0.027±0.011 indi- 
cates that the temperature fluctuations are ~16 % inside neb- 
ula. The temperature fluctuation in BoBn 1 is low, compared 
with typical Galactic PNe, i.e. f 2 <0.1 (cf. Zhang et al. 2004). 

When we take the temperature fluctuation effect into 
account, the derived ADFs become lower than the f 2 =0 
case; ADFs are 0.87±0.25 for C 2+ , 1.12±0.38 for C 3+ , 
2.95±1.32 for N 2+ , 0.32±0.16 for N 3+ , 2.63±0.55 for 2+ , 
and 1.69±0.37 for Ne 2+ , respectively. The great improve- 
ments of C 2+ and C 3+ might have been caused by the 
large temperature dependency of Cm] AA1906/09 and Civ 
AA1549/50; the energy difference between upper and lower 
level, AE = kAT, where AT = 75 380 K and 44820 K, re- 
spectively. This also implies that C 2+ and C 3+ abundances 
from ORLs are more reliable than those from CELs. Con- 
cerning 2+ ADFs, large discrepancy still exists. Since com- 
pared with the Cm] AA1906/09 and Civ AA1549/50 lines, 
the [O III] nebular lines depend more weakly on the electron 
temperature (A 7^14 100 K), so we realize that the temper- 
ature fluctuation model alone cannot improve CEL 2+ over 
> 1 dex. We need to seek other explanations for the large 
ADF(0 2+ ). Potentially this model can explain the discrep- 
ancies of N 2+ and Ne 2+ abundances, taking into account the 
uncertainties of measured fluxes of the observed ORLs N II 
and Ne II. 

In Table 17, we present elemental abundances from the 
CELs and ORLs, taking into account the above temperature 
fluctuations. The CEL C, N, and Ne abundances become com- 
parable to the ORL abundances. For the O abundance, there 
still exists a large discrepancy. 

4.1.2. High density components 

It was proposed by Rubin (1989) and Viegas & Clegg 
(1994). It assumes that small high density components within 
nebula weaken the intensity of the nebular lines due to their 



U7±0.01 
i.77±0.16 
i.66±0.21 
>.46±0.10 

i.35±0.14 



TABLE 16 

t 2 AND T FOR EACH ZONE. 

zone t 2 7q 

0.027±0.011 8920±840 

1,2 0.027±0.011 11540±400 

3 0.027±0.011 12 220±630 

4 0.027±0.011 12 950±610 

5 0.027±0.011 12 870±730 

6 0.027±0.011 12 790±840 

collisional de-excitation, assuming that chemical abundances 
are homogeneous. In this situation, the nebular to auroral 
line intensity ratios become smaller, from which we would 
derive falsely high electron temperatures. Accordingly, the 
CEL ionic abundances would be underestimated. Since the 
ORLs have very large critical densities, the ORL ionic abun- 
dances are hardly affected by collisional de-excitation. In the 
case of the halo PN DdDm 1, Otsuka et al. (2009) discussed 
the possibility that the 2+ abundance discrepancy could be 
explained by this model. However, as we argued in §3.3.2, 
no such high density components in BoBn 1 gas. Hence this 
model would not provide a sound ground for the observed 
abundance discrepancy. 

4.1.3. hydrogen deficient cold components 

This model was proposed by Jacoby & Ford (1983), Liu et 
al. (2000), Pequignot et al. (2002), Wesson et al. (2003), 
and others. This model assumes the situation as follows; the 
central star of a PN first ejects an envelope at low expansion 
velocity (of the order 10 km s" 1 ) with "normal" heavy metal 
abundances, and later ejects the high-velocity, hydrogen defi- 
cient, cold, and rich heavy metal components. Here, the ORLs 
are assumed to be emitted mainly from high-velocity hydro- 
gen deficient cold components, whereas the CELs are from 
the hot, normal metal gas surrounding the ORL emitters. 

So far, such components are directly or indirectly observed 
in the some PNe, for example, Abell 30 (Wesson et al. 2003), 
NGC 6153 and NGC 7009 (Barlow et al. 2006). In high- 
velocity components of Abell 30, Wesson et al. (2003) found 
that the electron temperature derived from O II lines is 500- 
2500 K and the ORL oxygen abundance is ~100 times larger 
than from CELs. Abell 30 is a well known PN, for having 
a hydrogen deficient central star and it is suspected to have 
experienced a very late thermal pulse. If this model is the 
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CELs ORLs CELs ORLs CELs ORLs 



TABLE 17 

The elemental abundances derived from CELs and ORLs in the case of t 2 yo. 



X/H 



log(X/H)+12" 



[X/H] 



He 




1.18(-1) ±2.12(-3) 




C 


1.44(-3) ± 3.31 (-4) 


1.44(-3) ±4.96(-4) 


9.16±0.10 


N 


1.77(-4) ± 5.37(-5) 


2.99(-4) ± 1.45(-5) 


8.25±0.14 


O 


6.35(-5) ± 7.38(-6) 


1.67(-4) ± 3.98(-5) 


7.80±0.05 


F 


5.43(-7) ± 1.57(-7) 




5.73±0.13 


Ne 


1.01(-4) ± 9.26(-6) 


1.71{-4) ± 4.08(-5) 


8.00±0.04 


S 


2.48(-7) ± 1.16(-7) 




5.32±0.22 


CI 


2.47(-9) ± 4.02(-10) 




3.39±0.07 


Ar 


2.36(-8) ± 3.61(-9) 




4.37±0.07 


Fe 


1.53(-7) ± 5.71(-8) 




5.18±0.17 


Kr 


<8.74(-10) 




<2.94 


Xe 


<1.33(-10) 




<2.12 


Ba 2.54(-10) ± l.OO(-lO) 




2.41 ±0.18 



11.07±0.01 
9.16±0.16 
8.48±0.23 
8.22±0.11 

8.23±0.11 



+0.17±0 
+0.77±0.11 +0.77±0 
+0.42±0.18 +0.66±0 
-0.89±0.07 -0.47±0 
+1.27±0.14 
+0.13±0.11 +0.36±0 
-1.80±0.23 
-1.94±0.09 
-2.18±0.10 
-2.29±0.17 

<-0.42 

<-().()9 
+0.23±0.18 



"The number density of the hydrogen is 12. 
b Solar abundances are taken from Lodders (2003). 

case of Abell 30, the ORL oxygen abundance might indicate 
the amount of O synthesized in the He-rich intershell. Barlow 
et al. (2006) measured the expansion velocities of [O III] and 
Oil lines in NGC 6153 and NGC 7009 and found that the 
Oil expansion velocity is smaller than [Olll]. Hence, they 
concluded that the O II and [O III] lines do not originate from 
material of identical physical properties. 

To verify whether the large O (and 2+ ) discrepancy in 
BoBn 1 is due to difference physical properties between O II 
and [O III] lines or not, following Barlow et al. (2006), we 
compare the expansion and the radial velocities of O II and 
those of [O III] lines. The resultant values are summarized in 
Table 18. The third and last columns are the radial velocity V r 
and twice the expansion velocity, 2V exp , respectively. We also 
estimated V r and 2V exp of Ne II and [Ne III]. In the last one of 
the line series of each ion, we present the adopted velocities 
with the bold face characters. These values are estimated from 
the line intensity weighted means. The radial velocities of the 
O II, [O III], Ne II and [Ne III] lines are almost consistent with 
the average radial velocity of 191.6 ± 1.3 km s" 1 that are de- 
rived from over 300 lines detected in the HDS spectra. How- 
ever, the 2y exp values of the O II and Ne II lines are^lOkm s" 1 
smaller than those of [O III] and [Ne III] lines, respectively. 
These findings does not agree to the foregoing high velocity 
hydrogen deficient cold model. We can attribute the differ- 
ence between the expansion velocities in ORLs and CELs to 
their thermal motions. Then, we can assume a presence of 
oxygen and neon-rich components with the same normal ra- 
dial velocity, surrounded by hot normal-oxygen and neon gas. 

Otsuka et al. (2008a) argued that the rich-neon abundance 
of BoBn 1 might have been caused by a late-thermal pulse, 
and the central star might have been hydrogen-deficient. At 
that phase, hydrogen-deficient, oxygen and neon-rich cold 
components might be incidentally ejected from the central 
star. The ORL overabundance could be an evidence indica- 
tive of relatively recent yields in the central star. Georgiev et 
al. (2008) found that the ORLs He, C, and O abundances in 
the nebula of NGC 6543 are in good agreement with those in 
the stellar wind zone, while Morisset & Gorgiev (2009) found 
that the ORL C, N, and O abundances in the nebula of IC 418 
are in good agreement with those in the stellar wind. In the 
halo PN K 648, Rauch et al. (2002) estimated the C, N, and O 
abundances using the stellar spectra; C=1.0(-3), N=1.0(-6), 
and O=1.0(-3). Their estimated C abundance agrees with the 



1.01 
1.16 
1.26 
U2 

M5 



TABLE 18 

Radial and twice the expansion velocities for 2+ and Ne : 



lines . 


Ion 






2Vexp 




(A) 


(km/s) 


(km/s) 


Oil 


4089.29 


+190.8 ± 2.8 


23.7 ± 6.4 




4153.30 


+206.8 ± 2.6 


33.6 ±5.6 




4292.21° 


+190.8 ± 2.8 


30.0 ± 4.9 




4638.86 


+190.1 ± 1.3 


31.0 ±2.2 




4641.81 


+187.2 ± 4.4 


24.2 ± 6.2 




4649.13 


+205.4 ± 0.9 


31.2 ±2.0 




4650.84 


+205.3 ± 1.6 


28.3 ± 1.7 




4673.73 


+185.0 ± 4.8 


55.2 ± 24.2 




4676.23 


+194.2 ± 5.8 


41.2 ± 17.3 




6721.39 


+194.2 ± 5.8 


24.5 ± 1.9 




Adopted 


+195.1 ± 3.0 


31.8 ± 6.9 


[OIII] 


4363.21 


+191.0 ±0.2 


40.3 ±0.1 




4931.80 


+194.2 ± 1.7 


47.5 ± 3.9 




4958.91 


+193.7 ±0.8 


40.5 ± 0.3 




5006.84 


+193.8 ± 1.3 


42.8 ± 0.3 




Adopted 


+193.7 ± 1.2 


42.2 ± 0.3 


Nell 


3694.21 


+193.8 ± 3.2 


34.4 ± 5.3 


[Ne III] 


3342.42 


+199.5 ±0.5 


42.7 ± 0.9 




3868.77 


+191.1 ±0.7 


41.5 ±0.1 




3967.46 


+191.9 ±0.1 


41.5 ±0.2 




Adopted 


+191.3 ± 0.6 


41.5 ± 0.1 



a O II A 4291.26,4291.86,4292.21,4292.98 are blended. 

ORL C abundance of 1.8(-3) within the errors (Otsuka 2007). 
To verify whether the ORL abundances in BoBn 1 indicate 
recent yields in the central star or not, we need to estimate the 
stellar abundances and compare those with the nebular ORL 
abundances using high-dispersion UV spectra. 

4.2. Comparison of elemental abundances from this work 
with others 

In Table 19, we compiled results for BoBn 1 from the past 
30 years. Our estimated CEL elemental abundances except 
for Fe are in good agreement with previous works. The large 
discrepancy for Fe between us and Kniazev et al. (2008) could 
be due to adopted electron temperatures for the Fe 2+ abun- 
dance estimation. The depletion of Fe relative to the Sun is 
almost consistent with that of Ar (Table 15). Since Ar and Fe 
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are not to be synthesized in low-mass stars, the abundances of 
these elements must be roughly same, if the large amount of 
the dust does not co-exist in the nebula. In addition, the mid- 
IR spectra show no astronomical silicate or iron dust such as 
FeS features. Therefore, our estimated Fe abundance seems 
more reliable than Kniazev et al. (2008). Sneden et al. (2000) 
estimated [Fe/H] = -2.37±0.02 as the metallicity of M 15 us- 
ing >30 giants. The [Fe/H] abundance of BoBn 1 corresponds 
to that of Ml 5 within error, implying that the progenitor of 
BoBn 1 might have formed at ~40 Gyr ago. 

BoBn 1 is a quite Ne-rich PN, and this Ne abundance is 
comparable with those of bulge and disk PNe; the averaged 
Ne abundance is 8.09 (CELs) and 9.0 (ORLs) for bulge PNe 
(Wang & Liu 2007) and 7.99 (CELs) and 9.06 for disk PNe 
(Tsamis et al. 2004; Liu et al. 2004; Wesson et al. 2005). 
The Ne isotope, 20 Ne is the most abundant, and it is not al- 
tered significantly by H- or He-burning (Karakas & Lattanzio 
2003). Therefore, the Ne overabundance would be due to an 
increase of the neon isotope, 22 Ne. During helium burning, 
14 N captures two a particles, and 22 Ne are produced. The Ne 
overabundance also implies that BoBn 1 might have experi- 
enced a very late thermal pulse (Otsuka et al. 2008a). If this 
is the case, Ne abundance might be an indirect evidence of 
He-rich intershell activity. The Ne overabundance would be 
concerned with extra mixing during the RGB phase, which 
would increase N. The Ne and N enhancements of BoBn 1 
might be also involved with the chemical environment where 
the progenitor formed. So far, 4 objects including BoBn 1 
have been regarded as Sagittarius dwarf galaxy PNe, and 3 
objects of them showed C, N, and Ne-rich ([C,N,Ne/0] > 0; 
cf. Zjilstra et al. 2006). 

The leading and trailing streams of the Sagittarius dwarf 
galaxy trace several globular clusters. Terzan 8 is a member of 
the Sagittarius dwarf galaxy. Mottini et al. (2008) investigated 
the chemical abundances in three red giants in Terzan 8. The 
averaged [Fe/H], [O/Fe], and [Mg,Si,Ca,Ti/Fe] among these 
objects were -2.37±0.04, +0.71 ±0.14, and +0.37±0.14, re- 
spectively. Note that the metallicity of Terzan 8 is very close 
to BoBn 1 . The amounts of O and other a elements such as 
Mg, Ne, S, and Ar do not significantly change during RGB 
phase. Therefore, the pattern of a elements derived from 
these RGB stars should be close to those in BoBn 1 's progen- 
itor. Based on this assumption, the initial [O/H] abundance 
of BoBn 1 is estimated to be -1.66±0.15, which is log(0/H) 
+ 12 = 7.03±0.15. [O/H]-+0.77 would need to have been 
synthesized in BoBn 1 during helium burning. Assuming that 
the 22 Ne(a,n) 25 Mg reaction is inefficient, the initial [Mg/H] 
abundance is estimated to be -2.07±0.15, which is compara- 
ble to the observed [S,Ar/H] in BoBn 1. The observed S and 
Ar abundances in BoBn 1 could show the original abundances 
of the progenitor. [Ne/H] = +2.16±0.18 could have been syn- 
thesized in BoBn 1 by 14 N capturing two a particles. Mottini 
et al. (2008) also estimated the light and heavy ^-process en- 
hancements; [Ba/Fe] = -0.09±0.17 and [Y/Fe] = -0.29±0.17. 
From these values, the initial ^-process elemental abundances 
would be -2.46±0.18 for heavy ^-process elements such as 
Xe and Ba and -2.66±0.18 for light s-process elements such 
as Kr. The progenitor of BoBn 1 could have synthesized [s/H] 
~+2 during He-burning phase. 

4.3. Comparison of observations and photo-ionization 
models 

To investigate the properties of the ionized gas, dust, and the 
PN central star in a self-consistent way, we have constructed a 



theoretical photo-ionization (P-I) model which aims to match 
the observed flux of emission lines and the spectral energy 
distribution (SED) between UV and mid-IR wavelength, us- 
ing Cloudy c08.00 (Ferland 2004). 

First, a rough value of the distance to BoBn 1 is necessary 
in fitting the observed fluxes. The distance to this object is es- 
timated to be in the range between 16.5 and 29 kpc (see Table 
1). Based on the assumption that BoBn 1 is a member of the 
Sagittarius dwarf galaxy, we fix the distance to be 24.8 kpc 
(Kunder & Chaboyer 2009). P-I model construction needs in- 
formation about the incident SED from the central star and 
the elemental abundances, geometry, density distribution, and 
size of the nebula. Once one gets a proper prediction of line 
intensities from a proper modeling procedure, the central stel- 
lar physical properties employed in the P-I model can give us 
a hint of the nebular evolutionary history or that of its progen- 
itor star. Especially, the central star temperature and the 
SED of the PN central star are an important physical parame- 
ter in constructing the correct P-I model. 

We estimated of 125 930 ± 6100 K using the energy bal- 
ance methods proposed by Gurzadyan (1997). Being guided 
by this T±, we used theoretical model atmosphere for a se- 
ries of values of T e ff to supply the SED from the central star. 
We used Thomas Rauch's non-LTE theoretical model atmo- 
spheres 7 for halo stars ([X,F] = and [Z] = -1) with the sur- 
face gravity logg = 6.0, 6. 125, 6.25, 6.375, 6.5, and 6.625. We 
varied T e ff and the luminosity to match the observations. 

For the elemental abundances X/H, we used the values from 
the case of t 2 = as a starting point. For the C, N, O, and Ne 
abundances, we used the CELs abundances. We assumed no 
high-density cold clumps. Using the HDS slit-viewer image 
(Fig. 1), we measured the radius of the outer nebular shell R out 
of ~0.6" (=0.072 pc) and fixed this value. We assumed the 
hydrogen density Nh to be a R~ 2 smooth distribution, i.e., Nn 
= Nu(Rm) x (Ri„/R) 2 . In the models, within a small range we 
varied X/H, Nu(Rm), and R{ n to match the observed line fluxes 
from UV to mid-infrared wavelength, including 2MASS JHK 
bands, and our mid-infrared bands, i.e., between 17 and 23 
/im (IRS B) and between 27 and 33 fjxa (IRS C). 

Since the IRS spectra show that the dust grains co-exist in 
the nebula of BoBn 1, we need information about the dust 
composition. Here we considered amorphous carbon and 
PAH grains. The optical constants were taken from Rouleau 
& Martin (1991) for amorphous carbon and from Desert et al. 
(1990), Schutte et al. (1993), Geballe (1989), and Bregman 
et al. (1989) for PAHs. The observed emission-line and base 
line continuum fluxes between 5.9 and 6.9 /mi (PAH 6.4 ^m) 
and between 7.4 and 8.4 /im (PAH 7.9 /xm) were used to de- 
termine the PAH abundance. We assumed that the gas and 
dust co-exist in the same sized ionized nebula. We adopted a 
standard MRN a" 3 5 distribution (Mathis, Rumpl & Nordsieck 
1977) with fl m i„=0.001 /im and a max =0.25 /im for amorphous 
carbon. For PAHs we adopted an a -4 size distribution with 
flmin=0. 00043 /im and a max =0.0011 /im. We adopted an R~ 2 
smooth dust density distribution. High-density clumped dust 
grains were not considered. 

In Table 20, we compare the predicted with observed rel- 
ative fluxes where /(H/3) = 100. For most CELs and He 
lines and the wide band fluxes, as well, the agreement be- 
tween the P-I model and the observation is within 30%. The 
poor fit of [Ni], [Oi], and [Sn] would be due to the as- 
sumed density profile. The relation between n e and LP (Fig. 

7 See http://astro.uni-tuebingen.de/~rauch/ 
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TABLE 19 

Elemental abundances derived by previous works and by this work. 



Abundances (log(X/H) + 12) 



Nebula Ret. 


He 


C 


N 


O 


F 


Ne 


S 


CI Ar Fe 


Kr 


Xe Ba 


BoBn 1 (1) 




9.02 


8.03 


7.74 5.85 7.96 


5.32 


3.39 4.33 5.08 


<2.88 


<1.97 2.30 


(2) 


11.07 


9.16 


8.49 


8.23 




8.22 










(3) 




9.16 


8.25 


7.80 


5.73 


8.00 


5.32 


3.39 4.37 5.18 


<2.94 


<2.12 2.41 


(4) 


11.07 


9.16 


8.48 


o.LL 




o.Zj 










(5)" 


11.11 


8.63 


7.96 


7.70 




1 Oft 

f.y\f 


5.01 


3.22 4.29 5.05 






(6)" 


1 1 .02 


9.20 


7.90 


7.70 




7.80 


5.80 


■ ■ ■ 4.70 • • • 






U) 


1 1 .00 


9.39* 


8.08 


8.03 




7.94 










(8f 


11.05 


8.95 


8.00 


7.83 




7.72 


5.50 


■■■ 4.50 ■•• 






(9) 


10.95 




7.70 


7.89 




8.10 


4.89 


... 4.19 ... 






(10) 














<5.45 








(11) 














5.67 








(12) 


10.98 


9.09 


8.34 


7.90 




8.00 




... 4.59 ... 






(13) 
















... 4.59 ... 






(14)" 


10.98 


8.48 


6.94 


7.70 




7.62 


6.48 








(15) 


11.06 




8.52 


7.89 




7.72 










(16) 


10.99 






7.44 




7.76 










(17) 


11.00 


9.20* 


7.64 


7.81 
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FIG. 12. — The predicted SED from the P-I modeling (black line) and the observed spectrum of BoBn 1 from the WE, Subaru/HDS, VLT/UVES, and 
SpitzerflRS (gray lines). The circles are 2MASS data. 



8 lower panel) shows that the electron density jumps from 
^2000 to ^6000 cm" 3 around [S II] emitting region. In our 
model, such a density jump was not considered. For important 
lines such as He I, II, and C III], [N II], [O II], [O III], [Ne III], 
fairly good agreements in calculating ionic abundances are 
achieved. However, most of the C and O ORLs fit to the ob- 



servations poorly. In most cases, the P-I models underesti- 
mate their line fluxes. As we discussed above, the O ORLs, 
and likely the Ne ORLs too, might be emitted from cold and 
metal-enhanced clumps. 

In Table 21, we list the derived parameters of the PN central 
star, ionized nebula, and dust, and in Fig. 12 we present the 
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TABLE 20 

Comparison between the PI-model and the observations. 



Ion/Band 


A kb 


type 


/(Cloudy) 


/(Obs.) 


Ion/Band 


Ai b 


type 


/(Cloudy) 


/(Obs.) 




(A/^im) 




[7(H/3)=100] 


[7(H/3)=100] 




(A/^tm) 




[/(H/3)=100] 


[Z(H/3)=100] 


He I 


4471 


ORL 


5.40 


4.82 


O II 


4094 


ORL 


1.35(-2) 


1 43(-21 


He I 


4922 


ORL 


1.33 


1.32 


On 


4152 


ORL 


7.18(-3) 


1.68(-2) 


He I 


5876 


ORL 


1.65(+1) 


1.81(+1) 


O II 


4294 


ORL 


6.05(-3) 


1.76(-2) 


He I 


6678 


ORL 


3.64 


4.02 


O II 


4651 


ORL 


5.61(-2) 


4.68(-2) 


Hell 


4686 


ORL 


2 79(411 


2.48(+l) 


[Ne II] 


12.81 


CEL 


9.33(-l) 


2.49 


rc ri 


8727 


CEL 


6.61(-2) 


7 931-21 


[Ne III] 


3343 


CEL 


1.42 


8 47(-l ) 


re tii 


2326 


CEL 


7.56(+l) 


3.51(+1) 


[Ne III] 


3869 


CEL 


2.74(42) 


2.17(42) 


Cm] 


1907 


CEL 


o 26(421 


8.39(+2) 


[Ne III] 


3968 


CEL 


8.25(41) 


6 .39(+l) 


Cm] 


1910 


CEL 


6.55(+2) 


6.03(+2) 


[Ne III] 


15.55 


CEL 


1.30(42) 


1.61(42) 


Civ] 


1548 


CEL 


3.75(+2) 


i 05(4-31 


[Ne III] 


36.01 


CEL 


1 lion 


1 33(+l I 

L.jj\-ri) 


Civ] 


1551 


CEL 


1.90(+2) 


5.19(+2) 


[Ne IV] 


2424 


CEL 


2.54(41) 


1 71(411 


Cll 


4267 


ORL 


3.34(-l) 


7 90C— 1 1 


[Ne IV] 


4725 


CEL 


L.JZ.y I) 


9.60(-2) 


ClI 


6580 


ORL 


4 84C-21 


3.75(-l) 


TNe vl 

[lit v j 


3346 


CEL 


1.65(-1) 


3 22(-11 


Cm 


4069 


ORL 


1 57C-1") 

L.J I \ L f 


2 49(-11 


TNe Vl 


3426 


CEL 


4 51C-11 


8 711-11 


C in 


4187 


ORL 


5.48(-2) 


3.46(-l) 


rs ni 


4070 


CEL 


5.2H-2) 


3.95(-l) 


Cm 


4649 


ORL 


1 48C-11 


2 6K-11 


rs ni 


4078 


CEL 


1.68(-2) 


2.45(-2) 


C in 


8197 


ORL 


5.60(-2) 


4 39(-l1 
°+.jy\ l) 


rs ni 


6716 


CEL 


1 01 (-11 

l.OH, LJ 


1 23(-11 


C IV 


4659 


ORL 


7.25(-3) 


1 19C-11 


rs ni 


6731 


CEL 


1.53(-1) 


2.16(-1) 


C iv 


7726 


ORL 


3 301-31 


3 1 1 (-21 


rs mi 
liij 


6312 


CEL 


9 57(-21 


4.80(-2) 


TNll 


5198 


CEL 


3.28(-2) 


2 73(-11 


rs mi 


9069 


CEL 




3.78(-l) 


TNll 


5200 


CEL 


1.26(-2) 


1 9K-11 


rs ml 


18.67 


CEL 


6.38(-l) 


6.92(-l) 




5755 


CEL 


1.38 


1.23 


rs ivi 

L<j i v j 


10.51 


CEL 


1.95 


1.92 


TN III 


6548 


CEL 


1.46(+1) 


1.56(+1) 


[CI III] 


5518 


CEL 


1 81 (-21 


1.80(-2) 




6584 


CEL 


4.31(+1) 


5.04(+l) 


rci ivi 

IK^l 1 V J 


8047 


CEL 


2.08(-2) 


2 I0(-21 


NlII] 


1750 


CEL 


7.43(+l) 


4.81(+1) 


[Ar III] 


5192 


CEL 


5 55(-31 


4.00(-3) 


N IV] 


1486 


CEL 


2.91(+1) 


4.58(+l) 


[Arm] 


7135 


CEL 


3.18(-1) 


2 73(-11 


Nil 


4176 


ORL 


3.26(-3) 


1.36(-2) 


[Ar III] 


7751 


CEL 


7.67(-2) 


6.101-2) 


Nil 


4239 


ORL 


1.56(-2) 


2.42(-2) 


I Ar ivi 


4711 


CEL 


6.40(-2) 


9.40(-2) 


Nil 


4435 


ORL 


9.53(-3) 


1.27(-2) 


[Ar IV] 


4740 


CEL 


6 82(-21 


9.00(-2) 


N II 


5005 


ORL 




1 7Q(_21 


[Ar IV] 


7171 


CEL 


2.34(-3) 


6 53(-31 


Nil 


5679 


ORL 


2.52(-2) 


1.62(-2) 


[Ar IV] 


7263 


CEL 


i 97(_31 

L.y I \ j) ) 


4.58(-3) 


N III 


41 10 


ORL 


1 0U-21 


1 40C— 21 


[Fe III] 


5271 


CEL 


3.89(-2) 


2 20l"-21 


NlII 


4379 


ORL 


4 44C-21 


5 97(-21 


[Fe III] 


4755 


CEL 


1 37(-21 

L.J 1 \ S.) 


2 IOf-21 


NlII 


4641 


ORL 


4 9ic^n 


3.46(-2) 


[Fein] 


4881 


CEL 


2. 13(— 2) 


2 10C— 2") 


roil 


5577 


CEL 


1.48(-3) 


1.70(-2) 


2MASS / 


1.24 




4.98(41) 


6.10(41) 


[Oi] 


6300 


CEL 


4.02(-2) 


8.72(-l) 


2MASS 7/ 


1.66 




3.09(41) 


5.50(41) 


[Ol] 


6363 


CEL 


1.28(-2) 


2.91(-1) 


2MASS £ 


2.16 




2.26(41) 


3.27(41) 


[On] 


3726 


CEL 


1.05(+1) 


1.09(+1) 


PAH 


6.40 




5.28(41) 


7.32(41) 


[On] 


3729 


CEL 


5.96 


6.61 


PAH 


7.90 




1.32(42) 


9.47(41) 


[On] 


7323 


CEL 


1.03 


1.02 


IRSB 


20.00 




6.08(41) 


6.48(41) 


[On] 


7332 


CEL 


8.21(-1) 


7.81(-1) 


IRSC 


30.00 




3.13(41) 


2.97(41) 


[Oiii] 


4363 


CEL 


6.86 


5.57 












[Om] 


4931 


CEL 


5.03(-2) 


4.36(-2) 












[Oiii] 


4959 


CEL 


1.23(+2) 


1.22(42) 












[OIII] 


5007 


CEL 


3.70(+2) 


3.51(42) 












[OIV] 


25.88 


CEL 


1.21(+1) 


1.25(41) 













TABLE 21 

The derived properties of the PN central star, ionized nebula, and dust by the P-I model. 



central star nebula dust 



d (kpc) 


24.8 


composition 


He:ll.ll,C:8.63,N:7.96,O:7.70, 


grains 


amC & PAHs 


U\L Q ) 


1180 


(conti.) 


F:5.85,Ne:7.90,S:5.01,Cl:3.22, 


M dust (M ) 


5.78(-6) 


T*(K) 


125 260 


(conti.) 


Ar:4.29,Fe:5.05,others:[X]=-2.13 


7dust (K) 


80-180 


log g (cm 2 s~ 


') 6.5 


/?in^out (") 


0.43/0.60 


Mdust/Mgas 


5.84(-5) 


composition 


[X,Y]=0, [Z]=-l 


N H (R in ) (cm- 3 ) 


3890 


Mdust («q yr" 


') ~3(-9) 


M corc (M ) 


~0.62 


geometry 


spherical 










logF(H/3) 


-12.44 










M gas /M atora (Mq) 0.09/0.04 
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FIG. 13.— The location of BoBn 1 and K 648 on the HR-diagram. The 
values of K 648 are from Rauch et al. (2002). The solid and broken lines 
represent the post-AGB He-burning evolution tracks of Vassiliadis & Wood 
(1994) for a metallicity of ~ 0.5 Z . 

predicted SED (black line). The predicted SED matches the 
UV to mid-infrared region well. Through the P-I modeling, 
we found the PN central star's parameters: r e ff = 125 260 ± 
200 K, U = 1 180 ± 240 L©, log g = 6.5, and a core mass of 
^0.62 Mq. The estimated 0.62 core mass is comparable to 
that of K 648 (0.62 Af©, Bianchi et al. 2001; 0.57 M© Rauch 
et al. 2002) and the high-excitation halo PN NGC 4361 (0.59 
M & : Traulsen et al. 2005). The ionized mass of BoBn 1 is 
0.09 M Q , which is comparable to K 648 (0.07 M©; Bianchi et 
al. 1995 ). In Fig. 13 we plot the locations of BoBn 1 and K 
648 (Rauch et al. 2002) and the post-AGB He-burning evolu- 
tionary tracks for LMC metallicity (Z^0.5Z©) by Vassiliadis 
& Wood (1994). These evolutionary tracks would suggest the 
possibility that the progenitors of BoBn 1 and K 648 were 
single 1-1.5 M© stars which would end their lives as white 
dwarfs with a core mass of ^0.6 M©. Alternatively, these 
halo PNe might have evolved from binaries composed of a 0.8 
Mq (= a typical halo star mass) secondary and a more mas- 
sive primary, and gained mass (^0.1 M©) from the primary 
through mass transfer or coalescence. 

For K 648, Alves et al. (2000) support a binary evolu- 
tion scenario. From F enhancement and similarity to carbon- 
enhanced metal poor stars (CEMP), Otsuka et al. (2008a) ar- 
gued that BoBn 1 might have evolved from a binary composed 
of a 0.8 Mq secondary and a > 2 M© primary star. We should 
consider two possibilities: these halo PNe have evolved from 
single stars or from binaries. 

The P-I model indicated elemental abundances except for C 
and S are in excellent agreement with those estimated by the 
empirical method using the ICFs. The discrepancy for C is 
due to the underprediction of C IV lines. The model predicted 
C + = 2.0(-5) and C 2+ = 3.4(-4), which are comparable to the 
observations. However, it predicted lower line fluxes of C IV 
AA1548/51 than the observations, and accordingly an under- 
estimated C 3+ as 6.95(-5). These might suggest that the ori- 
gin of C IV lines is not the ionized nebula but the stellar wind 
zone. The discrepancy for S could be due to the low ionic the 
S + abundance, which depends strongly on the assumed radial 
density profile. 

For BoBn 1, we have for the first time estimated a dust 
mass of 5.78(-6) M© and the temperature of 80-180 K. The 
dust in BoBn 1 is carbon rich. The dust composition suggests 
that BoBn 1 had experienced the TDU during the latest ther- 



mal pulsing AGB phase (TP-AGB). Since the TDU efficiently 
takes place in > 1 M© stars, the progenitor of BoBn 1 might 
be 1-3.5 Mq from the aspect of elemental abundances and 
dust composition. 

The dust-to-gas mass ratio ip of 5.84(-5) is much lower than 
the typical value in PNe (<-10 -3 ; Pottash 1984). For AGB 
stars, Lagadec et al. (2009) argued that t/j scales linearly with 
the metallicity. They assume the ratio is 

V> = ipQ x 10 [Fe/H] (7) 

where ip Q is 0.005. When we adopt [Fe/H] = -2.22 for BoBn 
1, we obtain -0 = 3.01(-5). 

Assuming that the inner and outer shells were expanding 
with ^10 km s" 1 and most of the dust was formed during the 
ending period of the thermal pulse AGB phase, we estimated 
the dust mass-loss rate Mj„ sl of ^3(-9) M© yr -1 . Lagadec 
et al. (2009) estimated Md ust in metal-poor ([Fe/H]~-1) car- 
bon stars IRAS 16339-03 17 and 18120+4530 in the Galactic 
halo and IRAS 12560+ 1656 in the Sgr stream. They estimated 
4-18(-9) M Q yr" 1 . For IRAS 12560+1656, Groenewegen et 
al. (1997) estimated M dus , of 1.9(-9) M© yr" 1 . These M du st 
values are comparable to that of BoBn 1 . 

4.4. Comparison of observations and theoretical models 

Through P-I modeling, we found two possibilities: BoBn 
1 might have evolved from (a) a 1-1.5 M© single star or (b) 
a binary composing of ^0.8 M© secondary and a more mas- 
sive primary. In this section, we explore these possibilities by 
comparing the observed and predicted elemental abundances 
employing theoretical nucleosynthesis models for low- to in- 
termediate mass stars. 

In Table 22, we present observed elemental abundances and 
predicted values from the theoretical models of Karakas & 
Lugaro (2009) for 1.0, 1.5, and 2.0 M© stars with Z=10" 4 
([Fe/H]^-2.3). The abundances from the models are the val- 
ues at the end of the AGB phase. For these models, Karakas & 
Lugaro (2009) chose an initial a-enhanced abundance pattern, 
i.e., [a/Fe]=+0.4. For ^-process elements, they chose scaled 
solar abundances, i.e., [X/Fe]=0. These [a/Fe] and [X/Fe] ra- 
tios are consistent with the values for RGB stars in Terzan 
8, therefore the assumption of initial abundances seems very 
reasonable for BoBn 1 . The accuracy of the predicted abun- 
dances by the models is within 0.3 dex. For the observed 
abundances, we adopted the 1 2 =0 CEL abundances except for 
C. The adopted C abundance was from ORLs. 

On the possibility (a) that BoBn 1 has evolved from a sin- 
gle star and has survived in the Galactic halo, the abundances 
of BoBn 1 except for N can be properly explained by the 1 .5 
Mq star model including a partial mixing zone of 0.004 M©, 
which produces a 13 C pocket during the interpulse period and 
releases free neutrons through 13 C(a,n) 16 0. This model as- 
sumes that the stars end as white dwarfs with a core mass of 
~0.7 Mq, which is comparable to our estimated core mass. 

The 13 C(a,n) le O reaction proceeds in the upper surface 
layer of the He-burning shell. The F and s-process ele- 
ments are synthesized by capturing these neutrons in the He- 
intershell. The observed F, probably Kr and Xe abundances 
are systematically larger (^+0.5 dex) than 1 .5 M© star + par- 
tial mixing model. This suggests that hydrogen mixing mass 
is likely to be > 0.004 M© or that BoBnl had experienced 
helium-flash driven deep mixing (He-FDDM, Fujimoto et al. 
2000; Suda et al. 2004) and obtained the extra neutrons. This 
process can occur in stars with [Fe/H] <-2.5 at the bottom 
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TABLE 22 

Comparison of observations and the theoretical models for 
single and binary stars withz=l(r 4 . 







Abundances (log(X/H) + 12) 




Model 


C 


N O F 


Ne Fe 


Kr 


Xe 


1.00 M 


8.13 


6.52 6.98 3.79 


6.45 5.18 


2.00 


0.97 


1.50 M Q 


9.21 


6.81 7.76 5.28 


8.11 5.21 


2.43 


1.46 


+ partial mixing 


9.17 6.78 7.85 5.30 


8.31 5.21 


2.47 


1.50 


2.00 M© 


9.55 


6.87 7.92 5.93 


8.66 5.24 


2.37 


1.53 


0.75 M + 1.50 M Q 


9.51 


7.09 7.67 5.04 


7.47 5.22 


2.11 


1.30 


0.75 M + 1.80 M 


9.31 


6.72 7.50 4.74 


7.17 5.20 


1.88 


1.07 


0.75 M +2.10M Q 


9.23 


6.58 7.41 4.69 


7.11 5.20 


1.73 


0.91 


BoBn 1 (? 2 =0) 


9.16 


8.03 7.74 5.85 


7.96 5.08 


<2.88 


<1.97 


K 648 (t 2 =0) 


9.25 


6.36 7.78 ... 


6.87 ... 







of the He-burning shell because the entropy barrier between 
the H- and He-shell becomes low. The lower limit to [Fe/H] 
for BoBnl is -2.46. This process would also produce 14 N 
through the 13 CQ?,7) 14 N reaction, by mixing protons into the 
He-burning shell. 

We note that BoBn 1 is similar to K 648, for the latter is also 
known as an extremely metal poor, C- and N-rich halo PN. On 
the assumption that K 648 has been evolved from a single star, 
the abundances of K 648 except for Ne can be explained by a 
1 .5 Mq model. K 648 would not have experienced He-FDDM 
for abnormally increasing N. 

However, can halo single stars with an initial mass of ~ 1 .5 
Mq and Z = 10 -4 survive up to now? Such stars would end as 
white dwarfs in ^2-3 Gyr. Their lifetime is much shorter than 
the age of Terzan 8; Forbes et al. (2004) estimated the age of 
Terzan 8 to be 13 ± 1.5 Gyr from an age-metallicity relation 
for the Sagittarius dwarf globular cluster. If the progenitor of 
BoBn 1 was a ^0.8 Mq single star, then it can have survived 
up to now, however it cannot evolve into a visible PN and can- 
not become extremely C-rich. To circumvent the evolutionary 
time scale problem, we should consider the other evolutionary 
scenario for BoBn 1 . 

4.5. The Origin of BoBn 1 

As mentioned earlier, BoBn 1 is likely to be a binary origin 
because the [C,N,F/Fe] abundances of BoBn 1 are comparable 
to those of the carbon-enhanced metal poor star (CEMP) HE 
1305+0132 (Schuler et al. 2007) and other CEMP stars. 

Most CEMP stars show large enhancements of C and N 
abundances. Some evolutionary models for CEMP stars have 
demonstrated that the C and N overabundances would be 
reproduced by binary interactions. Schuler et al. (2007) 
concluded that HE 1305+0132 might have experienced mass 
transfer and that s-process elements should be enhanced. Lu- 
garo et al. (2008) concluded that HE 1305+0132 consisted 
of ^2 M Q (primary) and ^0.8 M (secondary) stars with 
Z = 10~ 4 and that the enhanced C and F could be explained 
by binary mass transfer from the primary star via Roche lobe 
overflow and/or wind accretion. In Fig. 14, we present the 
diagram of [Xe,Ba/Ar]-[C/Ar]. The [Xe/Ar] for BoBn 1 is 
an upper limit. The data for Galactic PNe except BoBn 1 are 
taken from Sharpee et al. (2007). The data for s-process el- 
ements enhanced CEMP (CEMP-s) with [Fe/H] > -2.5 and 
C-rich AGB are from the SAGA data base (Suda et al. 2008). 
For PNe, we use Xe as a heavy s-process element and adopt 
Ar as a metallicity reference. For CEMP-.? and C-rich AGB 
stars, we use Ba and Fe as a metallicity reference. The di- 
agram indicates that C and s-process elements are certainly 



synthesized in the same layer and brought up to the stellar 
surface by the TDU. We note that the enhancement of heavy 
i-process elements in BoBnl is comparable to CEMP-.? stars, 
in particular CS22948-027 (Aoki et al. 2007;[Fe/H] = -2.21, 
[Ba/Fe] = +2.31, [C/Fe] = +2.12, [N/Fe] = +2.48). The chem- 
ical similarities between BoBn 1 and CEMP-s stars suggest 
that this PN shares a similar origin and evolutionary history. 

BoBn 1 is similar to K 648 in elemental abundances and 
nebular shape (see Fig. 1 and Table 19). K 648 has been for 
a long time suspected to have experienced binary evolution. 
Rauch et al. (2002) and Bianchi et al. (2001) analyzed the 
spectrum of the central star and estimated a core mass ^0.6 
Mq. The mass ^0.6 Mq corresponds to the initial mass of 
1-1.5 Mq from the HR-diagram as shown in Fig. 13. The ini- 
tial mass of 1-1.5 Mq suggests that K 648 might have evolved 
from a binary and accreted a part of the ejected mass by a 
massive primary or coalescence during the evolution. Alves 
et al. (2000) argued that K 648 has experienced mass augmen- 
tation in a close-binary merger and evolved as a higher mass 
star to become a PN. Such a high mass star would be a blue 
straggler. Ferraro et al. (2009) observed stars in the globular 
cluster M 30 using the HST/WFPC2 and concluded that blue 
stragglers are results of coalescence or binary mass-transfer. 

In view of the internal kinematics and chemical abun- 
dances, the progenitor of K 648 seems to be a binary. K 
648 has a bipolar outflow (Tajitsu & Otsuka 2006) and bipolar 
nebula (Alves et al. 2000). The statistical study of PN mor- 
phology shows that bipolar PNe have evolved from massive 
stars with initial mass > 2.4 Mq or binaries. Otsuka et al. 
(2008b) showed that the [C/Fe] and [N/Fe] abundances of K 
648 are compatible with CEMP stars. So far, there have been 
no reports on the detection for any binary signatures in both 
objects. The contradiction to the evolutionary time scale of 
this object can be avoided if BoBn 1 has indeed evolved from 
a binary. Similar to K 648, BoBn 1 could have evolved from 
a binary and undergone coalescence to become a PN. 

We explore the possibility of binary evolution of BoBn 1 
using binary nucleosynthesis models by Izzard et al. (2004, 
2009). We assume a binary system composed initially of 
a 0.75 Mq secondary and a 1.5/1.8/2.1 M primary with 
Z=10" 4 , separation = 219/340/468 R Q , respectively. We set 
eccentricity e = 0, common envelope efficiency a = 0.5 and 
structure parameter Ace = 0.5. He-FDDM is not considered. 
We set a 13 C pocket mass of 7.4 x 10" 4 Mq. The 13 C pocket 
contains 4.1 xlO" 6 Mq 13 C and 1.3xlO" y M 14 N. We set 
the 13 C pocket efficiency = 2 and adopt wind mass-loss rates 
from Reimers formula on the RGB and by Vassiliadis & Wood 
(1993) on the TP AGB. 

When we choose these initial primary and secondary 
masses and a, the binary system will experience Roche lobe 
overflow; and it will merge into a 1.2-1.4M single star at the 
TP- AGB phase and end its life as a white dwarf with a core 
mass of 0.62-0.68 Mq, 10.4-12.5 Gyr after the progenitor was 
born. We present the results in Table 22. The binary models 
might seem to explain the elemental abundances of BoBn 1 
and K 648. Concerning K 648, the 0.75 M & + 1 .5 M /l .8 M 
models fairly well match the prediction to the observed abun- 
dances. Among the models, the 0.75 Mq + 1.5 Mq model 
seems to be the best fit to BoBn 1 for the moment because 
this model can explain not only the abundance patterns but 
also the observed core mass (the predicted core mass ~0.64 
Mq). If the progenitor has experienced extra mixing in the 
RGB phase and increased N, the N overabundance can be ac- 
commodated by this model. The issues with the evolutionary 
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FIG. 14.— [Xe or Ba/Ar]-[C/Ar] diagram. The [Xe/Ar] value of BoBn 1 
is upper limit. For PNe, we adopt Ar as a metallicity reference. For CEMP-.V 
and C-rich AGB stars, we use Fe as a metallicity reference. 

time scale and C and N enhancements might be resolved si- 
multaneously if BoBn 1 has evolved from such a binary. At 
the present, we conclude consider that binary evolution sce- 
nario is more plausible for BoBn 1 . 

To discuss further the evolution of BoBn 1 and K 648, we 
need to increase detection cases of s-process elements and to 
investigate the isotope ratios of 12 C/ 13 C, 14 N/ 15 N, 16 0/ 17 0, 
and 16 0/ 18 0, which would be useful to investigate nucleosyn- 
thesis in the progenitors. It would be also necessary to com- 
pletely trace mass-loss history to improve mass-loss rate. So 
far, mass-loss history of evolved stars has been revealing by 
investigating spatial distribution of dust grains and molecular 
gas using far-infrared to millimeter wavelength data. The At- 
acama Large Millimeter Array (ALMA) and the thirty meter 
telescope (TMT) could open new windows to study the evo- 
lution of metal-poor stars such as halo PNe. 

5. CONCLUSION 

We have performed a comprehensive chemical abundance 
analysis of BoBn 1 using IUE archive data, Subaru/HDS 
spectra, VLT/UVES archive data, and Spitzer/IRS spectra. 
We calculated the ionic and elemental abundances of 13 el- 
ements using ORLs and CELs. The estimations of C, N, O, 
and Ne abundances from the ORLs and Kr, Xe, and Ba from 
the CELs are done for the first time. The C, N, O, and Ne 
ORL abundances are systematically larger than those from 
CELs. We investigated the cause of the abundance discrep- 
ancies. The discrepancies except for O could be explained by 
a temperature fluctuation model, and that of O might be due 
to hydrogen deficient cold components. 

In the optical high-dispersion spectra, we detected 
emission-lines of fluorine and s-process elements such as ru- 
bidium, krypton, xenon, and barium. The values of [F/H], 
[Kr/H], and [Xe/H] suggest that BoBn 1 is the most F-rich 



Acker, A., Marcout, J., Ochsenbein, F., Stenholm, B., & Tylenda, R. 1992, 
"Strasbourg - ESO catalogue of galactic planetary nebulae. Part 1; Part 
2" 

Aldrovandi, S. M. V. 1980, Ap&SS, 71, 393 

Aller, L. H. 1984 "Physics of Thermal Gaseous Nebulae", Astrophysics & 

Space Science Library vol. 112 
Alves, D. R., Bond, H. E., & Livio, M. 2000, AJ, 120, 2044 



among F detected PNe and is a heavy s-process element rich 
PN. The enhancement of C, N, and heavy ^-process is compa- 
rable to CEMP-s stars with [Fe/H] >- 2.5. This suggests that 
BoBn 1 shares a similar origin and evolutionary history with 
CEMP-s stars. 

We built photo-ionization model using non-LTE theoretical 
stellar atmosphere models to check consistency between el- 
emental abundances derived by empirical methods and from 
the model and to investigate the properties of the central star, 
ionized nebula, and dust in a self-consistent way to fit the 
IR wavelength region. In the modeling, we considered the 
presence of dust. We compared the observed elemental abun- 
dances with theoretical nucleosynthesis model predictions for 
single stars and binaries with Z = 10 -4 . The observed elemen- 
tal abundances except for N could be explained either by a 
1 .5 Mq single star model or a binary model composed of 0.75 
Mq + 1 .5 Mq stars. Through the modeling, we estimated the 
luminosity and effective temperature and surface gravity of 
the central star and the total mass of ionized gas and dust and 
even for the SED of BoBn 1 . Using theoretical evolutionary 
tracks for post- AGB stars, we found that the progenitor of the 
central star was perhaps a 1-1.5 Mq star and evolved into a 
system of a white dwarf with a core mass of ^0.62 Mq and 
an ^0.09 Mq ionized nebula. We estimated a total mass of 
5.8xl0" 6 Mq in the nebula, which composes of amorphous 
carbon and PAHs. The presence of carbon dust indicates that 
BoBn 1 has experienced the third dredge-up during the ther- 
mal pulse AGB phase. 

The progenitor might have been initially quite N-rich. The 
He-flash driven deep mixing might be responsible for the 
overabundance of N. From careful consideration of observa- 
tional results and a comparison between BoBn 1 and K 648 
in M 15, we propose that the progenitor was a 0.75 Mq + 1.5 
Mq binary with, e.g. an initial separation of 219 Rq and had 
experienced coalescence during its evolution to become a C- 
and N-rich PN. The similar evolutionary scenario would be 
also applicable to K 648. 
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TABLE A 1 

Observed and reddening corrected line ratios [7(H/3)=100] and identifications of BoBn 1. The few lines from the Table Al. 



A„bs Ion Ai a b Comp. /(A) /(A) 57(A) Source Note 
(A) (A) 



3301.48 


O III 


3299 


,39 


1 


0.423 


0.632 


0.031 


UVES1 


3314.43 


O III 


3312 


.33 


1 


0.418 


1.644 


0.038 


UVES1 


3336.76 


Ne ii 


3334 


,87 


1 


0.411 


0.104 


0.019 


UVES1 


4717.13 


[Ne IV] 


4714 


,25 


1 


0.041 


0.055 


0.003 


HDS 


4718.61 


[Ne IV] 


4715 


,80 


1 


0.041 


0.019 


0.002 


HDS 


4727.07 


[Ne IV] 


4724 


,15 


1 


0.039 


0.051 


0.002 


HDS 


4728.51 


[Ne IV] 


4725 


,62 


1 


0.038 


0.045 


0.003 


HDS 


4743.20 


[Ar IV] 


4740 


,17 


1 


0.034 


0.090 


0.004 


HDS 


4757.92 


[Fe III] 


4754 


,69 


1 


0.030 


0.021 


0.005 


HDS 


4788.04 


C IV? 


4785 


,90 


1 


0.021 


0.044 


0.005 


HDS 


4790.11 


N IV 


4786 


,92 
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0.020 


0.013 


0.005 


UVES1 


4792.87 


[F II] 


4789 


,45 
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0.020 


0.056 


0.005 


HDS 


4801.07 


[Fe II] 


4798 


,27 
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0.017 


0.008 


0.005 


HDS 


4805.62 


Ne II 


4802 


,58 
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0.016 


0.034 


0.006 


HDS 


4855.14 


He II 


4852 


,00 


1.00 


0.003 


0.007 


0.004 


HDS 


4856.93 


[Fe III] 


4853 


,70 


1 


0.002 


0.004 


0.003 


HDS 


4862.36 


He II 


4859 


,32 


1 


0.001 


0.869 


0.130 


HDS 


4864.46 


H4 


4861 


.33 


1 


0.000 


100.000 


0.202 


HDS 


4872.75 


[F II] 


4868 


,99 


1 


-0.002 


0.013 


0.003 


HDS 


4879.17 


[Co VI] 


4876 


,26 


1 


-0.004 


0.009 


0.004 


HDS 


4884.34 


[Fe III] 


4881 


,00 


1 


-0.005 


0.021 


0.005 


HDS 


4924.71 


He I 


4921 


,93 


1 


-0.016 


0.277 


0.011 


UVES1 


4925.21 


He I 


4921 


,93 


2 


-0.016 


1.043 


0.013 


UVES1 










Tot. 




1.320 


0.017 


UVES1 


4934.42 


[O III] 


4931 


.23 


1 


-0.019 


0.044 


0.004 


HDS 


4937.32 


Ba II 


4934 


.08 


1 


-0.019 


0.006 


0.002 


UVES1 


4961.87 


[O III] 


4958 


,91 


1 


-0.026 


29.603 


2.389 


HDS 


4962.03 


[O III] 


4958 


,91 


2 


-0.026 


32.311 


1.991 


HDS 


4962.28 


[O III] 


4958 


,91 


3 


-0.026 


60.490 


5.079 


HDS 










Tot. 




122.404 


5.956 


HDS 


4967.84 


C II 


4964 


.73 


1 


-0.027 


0.027 


0.004 


HDS 


5004.63 


N II 


5001 


,47 


1 


-0.036 


0.018 


0.003 


HDS 


5009.85 


[O III] 


5006 


,84 


1 


-0.038 


55.674 


13.550 


HDS 


5009.97 


[O in] 


5006 


,84 


2 


-0.038 


137.072 


3.874 


HDS 


5010.25 


[O III] 


5006 


,84 


3 


-0.038 


158.090 


16.602 


HDS 










Tot. 




350.836 


21.776 


HDS 


5018.53 


He I 


5015 


.68 


1 


-0.040 


0.291 


0.012 


HDS 


5019.01 


He I 


5015 


,68 


2 


-0.040 


1.481 


0.014 


HDS 










Tot. 




1.772 


0.019 


HDS 


5033.51 


[Fe IV] 


5030 


.33 


1 


-0.043 


0.089 


0.008 


HDS 


5035.14 


C II 


5032 


,07 


1 


-0.044 


0.053 


0.005 


HDS 


5038.93 


[Fe II] 


5035 


,48 


1 


-0.045 


0.011 


0.004 


HDS 


5039.16 


C II 


5035 


,94 


1 


-0.045 


0.018 


0.003 


UVES1 


5050.58 


He I 


5047 


.74 


1 


-0.048 


0.044 


0.009 


HDS 


5051.11 


He I 


5047 


,74 


2 


-0.048 


0.211 


0.011 


HDS 










Tot. 




0.254 


0.014 




5059.25 


Si II 


5055 


,98 


1 


-0.050 


0.007 


0.004 


UVES2 


5069.55 


Nil 


5066 


,46 


1 


-0.052 


0.042 


0.009 


UVES2 


5117.37 


O V 


5114 


,06 


1 


-0.063 


0.013 


0.030 


HDS 


5122.44 


[Fe III] 


5119 


,97 


1 


-0.065 


0.014 


0.004 


HDS 


5125.06 


C II 


5121 


,83 


1 


-0.065 


0.029 


0.002 


HDS 


5130.66 


Fe III 


5127 


,39 


1 


-0.066 


0.007 


0.002 


HDS 


9074.84 


[S HI] 


9068 


,60 


1 


-0.594 


0.378 


0.010 


UVES2 


9101.95 


Nil 


9096 


,16 


1 


-0.596 


0.011 


0.004 


UVES2 


9114.25 


Hell 


9108 


,54 


1 


-0.597 


0.032 


0.003 


UVES2 


9205.28 


Fe I 


9199 


,45 


1 


-0.603 


0.012 


0.002 


UVES2 


9216.31 


He I 


9210 


,34 


1 


-0.604 


0.093 


0.006 


UVES2 


9219.44 


He I 


9213 


,24 


1 


-0.604 


0.025 


0.004 


UVES2 


9223.29 


Mg II? 


9218 


,25 


1 


-0.604 


0.008 


0.003 


UVES2 


9230.87 


Hell 


9225 


,23 


1 


-0.605 


0.032 


0.006 


UVES2 


9234.92 


Hi 


9229 


,01 


1 


-0.605 


2.545 


0.043 


UVES2 


9350.69 


Hell 


9344 


,94 


1 


-0.613 


0.572 


0.011 


UVES2 



Raman line? 



